Preprint typosot in JHEP stylo - HYPER VERSION 



arXiv:1206.5225 [hep-th] 



Entanglement Entropy for Singular Surfaces 

^ 

O 
(N 

^ Robert C. Myers" and Ajay Singh''^ 

a Perimeter Institute for Theoretical Physics, Waterloo, Ontario N2L 2Y5, Canada 
OO * Department of Physics & Astronomy and Guelph- Waterloo Physics Institute, 

University of Waterloo, Waterloo, Ontario N2L 3G1, Canada 

Oh 

(D 

c"! Abstract: We study entanglement entropy for regions with a singular boundary in higher 

dimensions using the AdS / CFT correspondence and find that various singularities make new 
CN universal contributions. When the boundary CFT has an even spacetime dimension, we find 

ly-^ that the entanglement entropy of a conical surface contains a term quadratic in the logarithm 

CN of the UV cut-off. In four dimensions, the coefficient of this contribution is proportional to the 

central charge c. A conical singularity in an odd number of spacetime dimensions contributes 
\^ a term proportional to the logarithm of the UV cut-off. We also study the entanglement 

entropy for various boundary surfaces with extended singularities. In these cases, similar 
universal terms may appear depending on the dimension and curvature of the singular locus. 
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1. Introduction 



Entanglement entropy has emerged as an interesting theoretical quantity for studies of quan- 
tum matter. For example, it allows one to distinguish new topological phases or different 
critical points [1, 2, 3]. Entanglement entropy has also been considered in discussions of 
holographic descriptions of quantum gravity, in particular, for the AdS/CFT correspondence 
[4, 5, 6]. In this context, as well as characterizing new properties of holographic field theories 
[7], it has been suggested that entanglement entropy may play a fundamental role in the 
quantum structure of spacetime, e.g., [8, 9]. 

In quantum field theory (QFT), the typical calculation of entanglement entropy begins by 
choosing a particular spatial region V and integrating out the degrees of freedom in the com- 
plement V . Then with the resulting density matrix pv, one calculates: See = — Tr (pv log pv)- 
Unfortunately, this calculation generically yields a UV divergent answer because the result is 
dominated short distance correlations in the vicinity of the boundary S = dV . However, if we 
regulate the calculation by introducing a short distance cut-off 5, the entanglement entropy 
exhibits an interesting geometric structure [10, 11]. For example, with a QFT in d spacetime 
dimensions, this allows us to organize the results as follows: 

■^EE = + H , (1-1) 

where each of the coefficients Ci involves an integration over the boundary S. In particular, 
the first two coefficients may be written as 



C2= d'^-'^yVhh2 = h2Aj:, (1.2) 

C4 = / d'^-^yVh\hi,in + hi,2K''Kl^] . (1.3) 

Of course, the leading term then yields the famous 'area law' result [12, 13]. The geometry of 
the boundary becomes even more evident in the second coefficient C4 with the appearance of TZ 
and K^^, the intrinsic Ricci scalar and the extrinsic curvature of this surface.^ The coefficients 
hk^a above will depend on the detailed structure of the underlying QFT. In particular, we 
note that they may become functions of the cut-off through their dependence on mass scales 
/ij in the QFT since the latter will only appear in the dimensionless combination fj,i6, i.e., 
hk,a = hk^a{fJ-i^) [11]- The geometric character of the entanglement entropy illustrated here 
naturally follows from the implicit choice of a covariant regulator in the QFT and the fact 
that the UV divergences are local. 

Unfortunately, the coefficients appearing in the expansion above are scheme dependent. 
Clearly, if we shift 6 — )■ a6, we find h^^a — ^ ^fc,a = tt^~'^^fc,a(oA'j'5)- Hence the regulator 
dependence here comes both from the implicit dependence on mass scales in the QFT and 
the 'classical' engineering dimension of the individual coefficients. Of course, the latter can 



'^We are assuming here that the background geometry is simply d-dimensional flat space. Otherwise addi- 
tional contributions could appear in C4 involving the background curvature. 
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be avoided by carrying the expansion in eq. (1.1) to sufficiently high orders. In particular, 
one may find a logarithmic contribution to the entanglement entropy 



= Cd \og{5/L) (1.4) 

where L is some (macroscopic) scale in the geometry of S. At this order in the expansion, 
the coefficient Cd is dimensionless. Further it may be natural to eliminate any intrinsic scales 
in the QFT by focussing on fixed point theories of the RG flow.^ In this case, shifting 5 as 
above makes a finite shift in the entanglement entropy but Cd is left unchanged. Hence one 
expects that this coefficient contains universal information that characterizes the underlying 
CFT. For example, in four dimensions, this universal coefficient is simply given by [16]^ 



Cd = — / (PyVh 
27r is 



(1.5) 



where a and c are the two central charges of the CFT. In principle then, this provides an 
approach to determine these central charges. For example, calculating the entanglement 
entropy for a sphere yields Cd oc o, while only c appears for a cylinder [16]. 

The preceding discussion implicitly assumes that the boundary of the region in question 
is smooth. The purpose of this paper is to explore the effects of geometric singularities 
in this 'entangling surface' S. In part, our motivation comes from the observation that in 
three dimensions, if the boundary contains corners or kinks, the corresponding entanglement 
entropy will contain additional logarithmic contributions [17, 18, 19, 20] 

= qsi^e) \og{5/L) , (1.6) 

where IS.9 is the opening angle of the kink - see figure 1. As a function of A0, the coefficient 
53 (A0) satisfies certain simple properties [17, 18, 19]. In particular, (73 (A0 = vr) = since 
the E becomes a smooth surface when the angle is set to vr. Strong subadditivity can be 
used to argue that in general g3(A0) must satisfy certain inequalities, e.g., q-i{/S.9) > and 
(^^{M) < [18, 19]. If the QFT is in a pure state, we have 5ee(^) = S'ee(V') and so 
q3{A9) = (73(271 — AO) in this case. Further, in examples [17, 18, 20, 19], one finds for a 
small opening angle: qsiAO — t- 0) oc 1/A^. However, we must add that the precise universal 
information contained in q^^AO) remains to be understood. 

A natural question is to ask whether similar contributions arise for singular entangling 
surfaces in higher dimensions. If yes, we can ask what the geometric dependence of these new 
terms is. Further, if we focus on CFT's, we might ask if the coefficients in these contributions 
have a simple dependence on the central charges, analogous to that in eq. (1.5). The AdS/CFT 
correspondence [21] provides a simple framework with which we may begin to address these 
questions. In fact, using the standard calculation of holographic entanglement entropy [4, 5], 



■^Universal terms may also appear in other circumstances, either as a finite contribution, e.g., [1, 3] or even 
when exphcit mass scales are present [14, 15]. 

^ Again, we are assuming that the background geometry is flat. 
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the logarithmic contribution (1.6) associated with a kink in three dimensions has already been 
identified in [19]. More generally, this approach allows us to easily study boundary CFT's 
in a variety of dimensions and further the geometric structure of the entanglement entropy 
becomes readily evident in holographic calculations [15, 22]. In the following then, we use 
holography to study some simple singular entangling surfaces in higher dimensions and we 
find a variety of new universal contributions. While these are just first steps towards a full 
understanding of these universal terms, our results indicate a rich geometric structure for the 
entanglement entropy of singular surfaces. 

The remainder of this paper is organized as follows: In the next section, we give a brief 
overview of our calculations and summarize our main results. In section 3, we consider 
entangling surfaces with a conical singularity for boundary CFT's with d = 4, 5 and 6. In 
these cases, the singularity in the geometry of the entangling surface is confined to a single 
point and so we broaden our calculations to consider extended singularities in section 4. 
There we find that the appearance of universal terms in the entanglement entropy depends 
on the dimension and the curvature of the singular locus. Section 5 presents calculations 
of holographic entanglement entropy for singular surfaces in boundary CFTs which are dual 
to the Gauss-Bonnet gravity. These calculations allow us to examine the dependence of the 
universal terms on the central charges of the underlying CFT. In section 6, we briefly discuss 
our results and consider future directions. Appendix A describes an alternate calculation 
of the entanglement entropy associated with a conical singularity. In particular, we make a 
conformal transformation from R'^ to Rx S'^~^ and so the conical entangling surface becomes 
a cylinder in the latter background. In appendix B, we give certain details for lengthy 
calculations presented in sections 3, 4 and 5. 

While we were in the final stages of preparing this paper, ref. [23] appeared which contains 
some results which overlap with ours. Particularly, authors have calculated the holographic 
entanglement entropy for a cone and a crease for a four-dimensional CFT. We also learned of 
an upcoming paper [24] where the entanglement entropy for a cone in a six-dimensional CFT 
is studied. 

2. Singular entangling surfaces and summary of results 

In the sections 3, 4 and 5, we will describe in detail various holographic calculations of 
the entanglement entropy for certain singular surfaces. Each of these calculations is quite 
lengthy and individually they are not very enlightening. Hence in this section, we provide 
an overview of these calculations and a summary of our results. We begin by describing the 
kinds of singular entangling surfaces which we will consider. 

Let us go back to eq. (1.6) for three dimensions. In this case, the entangling surface is a 
one-dimensional curve and the 'singular surface' would be one containing a kink or a corner 
where the direction of the tangent vector changes discontinuously at a point. We can charac- 
terize this behaviour by saying that the geodesic curvature of the curve contains a (5-function 
singularity. In higher dimensions, the entangling surface is a (d — 2)-dimensional submanifold 
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(a) 



(b) 



(c) 



Figure 1: (Colour Online) Panel (a) shows a kink in constant Euclidean time slice in d = 3. Panel 
(b) shows the cone ci. Panel (c) shows the crease k x R^. 

which divides into two a time slice of the d-dimensional background spacetime. In this case, 
the natural extension of geodesic curvature is the extrinsic curvature of these surfaces. How- 
ever, a distinct feature characterizing the geometry of these higher dimensional surfaces is 
now their intrinsic curvature. Of course, for a fixed background, these two curvatures will be 
related to each other (and the background curvature) through the Gauss-Codazzi equations. 
However, it is worth noting that in discussing singular surfaces, we might consider singular- 
ities in either the extrinsic curvature or the intrinsic curvature. In particular, as we show 
with simple examples below, it is possible to construct surfaces where the intrinsic curvature 
is everywhere smooth while the extrinsic curvature is singular. The other possibility which 
we consider is when both the extrinsic and intrinsic curvatures have singularities.^ In either 
case, the examples which we consider below contain (5-function singularities. That is, in all 
of our examples, the curvatures characterizing the entangling surface are finite and smooth 
everywhere, except for a particular locus or subset of points. 

In considering such singular entangling surfaces, we introduce an intuitive nomenclature 
to simplify the discussion: kink, cone and crease. Examples of each are illustrated in figure 
1. To explain these terms, it suffices to consider evaluating the entanglement entropy in flat 
d-dimensional background R'^. With Euclidean signature, the metric can then be written in 
'cylindrical' coordinates as 



*We will assume the geometry of the background is everywhere smooth in the following. Combined with 
the Gauss-Codazzi equations, this assumption rules out the possibility that the extrinsic curvature may be 
smooth while the intrinsic curvature is singular. 



m 
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where denotes the round metric on a unit n-sphere. Hence we have d = 3 + n + m. The 
entanghng surface will be a (d — 2)-dimensional geometry embedded in this background (on 
a surface of constant t^). 

The natural model of a kink is then given in d = 3 (n = = m) as the union of two rays: 
k = {^E = 0, p = [0, oo), ^ = ibQ}. With this choice of angles, the opening angle between the 
two rays is = 2Q. Similarly, our prototype for a cone is given in d = 3 + n as the surface: 
Cn = {^E = 0, p = [0, oo), 9 = n}. Hence with this construction, we are only considering cones 
with a spherical cross-section S*". We define a crease to be the higher dimensional extension 
of either of these singular surfaces where we take the direct product of a kink A; or a cone c„ 
with some other manifold. With the flat background (2.1) above, the natural extension of 
the previous constructions gives the crease k x i?"* in d = 3 + m or a conical crease Cn x i?™ 
in d = 3 + n + m. 

In a certain sense the kink is simply a lower dimensional version of a cone, i.e., k = 
cq. However, we prefer to distinguish these two classes to highlight the difference in the 
singularities noted above. If we consider the crease k x i?*", the intrinsic geometry of this 
submanifold is smooth everywhere including p = and in fact for the construction above, the 
intrinsic curvature vanishes everywhere. On the other hand, the extrinsic curvature contains 
a (5-function singularity at the tip of the kink, i.e., p = 0. Again for the above construction, 
the extrinsic curvature vanishes everywhere away from this singularity. Now considering the 
conical crease c„ x i?™, one finds that both the intrinsic and extrinsic curvatures have 6- 
function singularities at the tip of the cone, i.e., p = 0. For the prototype constructed above, 
the intrinsic curvature vanishes away from this singular point but the extrinsic curvature is 
nonvanishing. 

In our calculations of the entanglement entropy, it will also be useful to consider the CFT 
in a background of the form i?^"^" x S*™ , for which we write the metric as 

ds^ = dtl + dp^ + p2 ((i02 ^ gjj^2 Q^^2^^ ^ ^2 ^^2^ _ (2.2) 

This allows us to consider creases of the form k x or c„ x S"". With these geometries, the 
intrinsic curvature of the entangling surface is nonvanishing everywhere but in particular, the 
geometry of the singular locus is now 5*". 

Now after this description of our singular surfaces, we briefly explain our holographic cal- 
culations and the results. We use the AdS/CFT correspondence to calculate the holographic 
entanglement entropy for a boundary CFT with the above singular entangling surfaces. In 
different calculations, we use the following bulk geometries: 

r2 / m \ 

ds^ = —\ dz^ + dtl + dp^ + {d9^ + sin^ ^d^l) + ^ [dx'f , (2.3) 

ds^ = ^ (dz^ + /i(z) [dtl + dp^ + p2 [de^ + sin2 OdQl)] + f2{z) . (2.4) 

Here L is the AdS curvature scale and z is the radial coordinate with the asymptotic bound- 
ary at z = 0. Of course, eq. (2.3) is just AdS^+i where d = 3 + n + m as above and with a 
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flat boundary metric written in the coordinates given in eq. (2.1). The geometry in eq. (2.4) 
approaches AdS^+i asymptotically with the curved boundary metric given in eq. (2.2). How- 
ever, because this boundary geometry is not conformally flat, the bulk geometry is not purely 
AdSrf+i and the functions, /i and /2, must be determined by evaluating the gravitational 
equations of motion. With these bulk metrics, (2.3) and (2.4), we can readily calculate the 
holographic entanglement entropy for the kink, cone or crease geometries described above. 

Recall the proposal for calculating holographic entanglement entropy [4, 5]: Having se- 
lected a particular entangling surface S in the d-dimensional boundary theory, we must find 
the minimal area surface m (of dimension d — 1) which extends in the bulk geometry such 
that the boundary of m matches S at the asymptotic AdS^+i boundary. The entanglement 
entropy is then given by^ 

See = -aZi min [^('^)] • (2-5) 

This conjecture passes a variety of consistency tests, e.g., see [5, 6, 25, 22] and in fact, for 
spherical entangling surfaces, there is a derivation of holographic entanglement entropy which 
yields this result (2.5) [26]. An implicit assumption in eq. (2.5) is that the bulk physics is 
described by (classical) Einstein gravity. However, this proposal was extended [22, 27] to 
certain dual gravity theories with higher derivative interactions, namely. Lovelock theories 
[28]. In this extension, which we use in section 5, one finds the surface which minimizes a new 
functional of the geometry of m. In particular, the latter was chosen to match an expression 
for the black hole entropy in the Lovelock theories [29]. 

The details of our holographic calculations of the entanglement entropy for kinks, cones 
and creases are given in the following sections. As expected, these calculations contain a 
variety of divergences that must be regulated. However, we pay particular attention to the 
question of whether there are new singular contributions associated with geometric singularity 
in the entangling surface. Further, in certain cases, we find that the singularity produces new 
universal contributions. Table 1 summarizes our results. 

As described above, if the entangling surface has a kink in d = 3, there is a new universal 
contribution (1.6) to the entanglement entropy. We review the holographic calculation of 
[19] which reveals this result in section 3. These calculations are extended to creases of the 
form k X R"^ in section 4 and there we find no universal contribution, i.e., no log 5 term, for 
these cases. However, we note that in this case the locus of the singularity, i.e., the tip of the 
crease, is flat. The results become more interesting if this locus is curved. If the singular locus 
is curved but odd dimensional, we again find no logarithmic contribution. However, if the 
singular locus is both curved and even dimensional, we find that the entanglement entropy of 
the crease contains a log 5 contribution. While we have a limited number of explicit examples 
of this behaviour, i.e., the case of A; x S^, it suggests to us that generally for creases of the 

^AU of our calculations use Euclidean signature for which the minimal area surface is appropriate. If the 
calculation is done in a Minkowski signature background, one seeks to extremize the area but this surface will 
only be a saddle point. In either case, the area must be suitably regulated to produce a finite answer. 
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d 


Background 


Geometry of 


Crease 


Crease 


Expected 


New 




spacetime 


entangling surface 


dimension 


curvature 


Divergences 


Divergences 


3 




k 





- 


1/5 


log(<5)* 


4 




Cl 





- 


1/52 , iog(5) 


WW 


5 




C2 





- 


1/5^, 1/5 


log(5)* 


6 




C3 





- 


l/5\ l/5\ log(5) 


log'(<5)* 


>3 


R'^ 


/C X i?'^"^ 


d-3 


flat 


l/5'^-2 


1/5^-^ 


4 


R^ X 51 


A: X 51 


1 


flat 


1/52 


1/5 


5 


R^ X 52 


A: X 52 


2 


curved 


l/5\ 1/5 


1/52, log(5)* 


6 


i?3 X 


A: X 53 


3 


curved 


l/5\ 1/52 


l/<53, 1/5 


6 


R^xS^ 


kx{R^ X 52) 


3 


curved 


1/52 


l/5\ 1/5 


5 




Cl X i?i 


1 


flat 


l/5\ 1/5 


\og{5)/5 


6 




Cl X 


2 


flat 


1/52, iog(5)* 


\og{5)/5 


7 




Cl X R^ 


3 


flat 


1/5^ , 1/(53, 1/5 


log(5)/53 


5 


R^ X 51 


Cl X 51 


1 


flat 




\og{5)/5 


6 


i?4x52 


Cl X 52 


2 


curved 


l/5\ 1/52, log 5 


\oi{5Y 


7 


i?4x53 


Cl X 53 


3 


curved 


1/53, 1/5 


log(5)/53 


6 


R^ 


C2 X i?^ 


1 


flat 


l/5\ 1/52, log(5)* 


1/5 


7 


R^ 


C2 X i?2 


2 


flat 


l/5^ 1/53, 1/5 


1/52 


6 


R^ X 51 


C2 X 51 


1 


flat 


1/54, 1/52^ i^g(^)* 


1/5 


7 


R^ X 52 


C2 X 52 


2 


curved 


1/55, 1/53, 1/5 


1/52, log(<5)* 



Table 1: Summary of the divergent terms for various singular surfaces from our holographic calcula- 
tions. Here d is the spacetime dimension of the CFT background, which can be both flat or curved. 
The 'expected' divergences are those which can arise with a smooth entangling surface - see discussion 
in the introduction. The 'new' divergences are produced by the singularity in the surface and vanish 
when the surface is smooth, i.e., = 7r/2. Any universal terms are marked with a - see the 
discussion in the main text. 



form k X Yj2m.^ there are new universal terms of the form 

/ d2™2/^/)^[7^^1og<5 (2.6) 

J a 

where a is the singular locus on the entangling surface and \RJ^\ represents some curvature 
invariant containing m powers of the curvature on this submanifold. 

We also consider entangling surfaces with conical singularities in section 3. In this case, 
if the boundary CFT lives in an odd number of spacetime dimensions, i.e., Cn with n even, 
we find that the singularity contributes a log 5 term to the entanglement entropy. However, 
for an even dimensional boundary theory, i.e., Cn with n odd, we find that the new universal 
contribution actually diverges as log2j. However, we note, and explain in detail in section 6, 
that part of this log25 contribution can be identified with the 'smooth' contribution given in 
eq. (1.5). That is, part of this divergence should be associated with correlations away from the 
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singularity and so depends on details of the smooth part of the geometry away from the tip of 
the cone. However, we also argue part of the contribution is associated with the singularity 
itself and so should still be expected to arise for more general situations, independent of 
this smooth geometry. Given that our boundary field theory is a CFT, we might ask if the 
coefficients of these new universal contributions are simple functions of the central charges. 
As a step in this direction, we work with Gauss-Bonnet gravity in the bulk in section 5, as 
this allows us to begin to distinguish the boundary central charges, e.g., see [30, 31]. In the 
case of even dimensions, we see that the coefficient of the log^S contribution is proportional 
to a particular central charge, i.e., for d = 4, it is c. However, for odd dimensions, the log (5 
term does not yield any such simple result. 

The holographic calculations are also extended to consider conical creases of the form 
Cn X i?*" or c„ X S*™ in section 4. For these cases, we find that the nature of universal 
contributions again depends on the dimension of both the full spacetime and the singular 
locus, as well as the curvature of the latter. In particular, if both of these dimensions is even 
and the locus is curved, e.g., ci x S^, then a log^S term arises. Alternatively, if the background 
is odd dimensional but the locus is even dimensional and curved, e.g., C2 x 5"^, then a log (5 
contribution appears. In any other cases, the singularity does not contribute any universal 
terms of this form to the entanglement entropy. Our results again suggest the appearance 
of universal contributions of the form given in eq. (2.6) for odd dimensional theories, while 
similar terms with a log^(5 divergence seem to be present for even dimensional CFT's. 

In discussing these results, we need to be careful about an important point. As we have 
illustrated with eq. (1.5), in even dimensional CFT's, the smooth part of the entangling 
surface will already produce a universal term proportional to log (5. Hence, as noted in our 
description of the results for cones, we must distinguish this term from universal contributions 
associated with the singularity. Other cases where this issue arises include: ci x i?^, C2 x 
and C2 X S^. All of these examples are in six dimensions where we do not have the analog 
of the d = 4 expression in eq. (1.5). So while a detailed comparison is not possible, in 
considering the corresponding holographic calculations in detail, we see that the coefficient 
of the log 5 receives contributions at all values of the radius p and that there is no singularity 
at p = 0. Hence we can clearly infer that this contribution is coming from the smooth part of 
the geometry and the singularity is not making a universal contribution to the entanglement 
entropy. 

To close this section, let us note that for many of the examples in table 1, there were 
no logarithmic terms in the entanglement entropy. In some of those cases, it may still be 
that the finite contribution exhibits some universal behaviour but we did not investigate this 
possibility here. 

3. EE for singular embeddings 

In this section, we will study entanglement entropy (EE) with singular entangling surfaces in 
a flat background for holographic CFTs which are dual to Einstein gravity. The simplest case 
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in this category is the kink in d = 3. It is already known that EE for a kink has a logarithmic 
divergence [17, 18, 20]. This calculation for holographic EE was first done by Hirata and 
Takayanagi in [19]. So before calculating EE for cones in higher dimensions, we briefly review 
this case. 

We begin with {d,n,m} = {3,0,0} in metric (2.3). The kink in the boundary is defined 
hy p £ [0,H] and 9 G where H is an IR cut-off. The holographic entanglement 

entropy for this geometry is given by (2.5), that is the area of the minimal area surface which 
hangs in the bulk and is homologous to the kink on the boundary. We assume that the induced 
coordinate for the minimal area surface are {p, 9) and the radial coordinate z = z{p, 9). Now 
we can find the induced metric h^^, over the surface and the entanglement entropy is given 
by 

= dpd9Vh = ^ j dpd9^^/p^ + p^z'^ + i\ (3.1) 

where z' = dpZ and z = dgz. Here we point out that the EE for entangling surface S in 
d-dimensional field theory will be represented by S'rfl^. Now we can easily find the equation 
of motion for z{p,9). From scaling symmetry of the AdS space and the fact that there is no 
other scale in the problem, we can say that 

z = ph{9) , (3.2) 

where h(9) is such that ^ — t- as — t- ±17. After using this ansatz, the entropy functional 
reduces to 

47rL2 dp f^~^ Jl + h^ + h? 

d9 , (3.3) 



^ J 5 /ho P Jo 

where we have introduced the UV cut-off at z = 6, h = dh/d9 and defined ho such that at 
^ = 0, h{0) = ho- Note that /i(0) = and ho is the maximum value of h{9). Also, e is a 
function of p defined using (3.2), such that at z = 5, h{0, — e) = 6 /p. Further, the substitution 
of ansatz (3.2) in equation of motion for z{p, 9) gives 

h{l + h^)h + 2h^ + {h^ + l)(/i2 + 2) = . (3.4) 

For this equation of motion, we can easily see that there exists a quantity K which is conserved 
along 9 translation and is given by 

. (3.5) 

where we have used /i(0) = and /i(0) = ho to get the expression on the right hand side. Now 
plan is to convert 9 integral in (3.3) to integral over h and then separate the divergent part 
in the integral. We also make the coordinate transformation to y = y^l/Zi^ — l/Zigi where 
y — )• oo as we approach the boundary. After this coordinate transformation, the integrand 
has following divergence in the limit y — t- oo: 



'l + hl{l + y^) ^ , , ^ 1 



2 + /l2(l + y2) 



~ 1 + O ( :3 ) . (3.6) 
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So now we can write EE as 

47rL2 '-H 



<5//io P Jo 




(3.7) 



which can be simphfied to give 



^3 



H ( b 

j+q3{n)log(- 



+ 0(6) 



where q^iO,) is 



dy 



1 



'l + /i2(l + y2) 



(3.8) 



(3.9) 



2 + hl{l + y^)_ 

Note that ^3(^1) is the cut-off independent term in the EE for the kink. After this quick review 
of the calculations by Hirata and Takaynagi, we turn towards the cone in higher dimensions 
where we will see log'^iJ divergence for even d. 

3.1 Cone in d = 4, 5 and 6 CFT 

In this section, we will calculate EE for cone Cn in some even and odd dimensional spacetime. 
We will give detailed calculations for ci in d = 4 and discuss the final results for C2 and C3 in 
d = 5, 6 dimensional CFTs^. However, in the beginning of the calculations, we will keep the 
discussion general for arbitrary d. 

With m = in the bulk metric (2.3), we can define the cone geometry by ^ G [0,^2] 
and p £ [0,H], where H is the IR cut-off for the geometry. Now, we define the minimal 
area surface, that gives the entanglement entropy, by coordinates (p, 0,^*) where .^"s are 
coordinates on sphere 5" in (2.3). As for the cone Cn, we have a rotational symmetry SO{n+l) 
along the sphere 5", the radial coordinate z will depend only on 6, i.e., z = z{p,6). Then, 
the induced metric is given by 



[1 + Z-) 



Z 



-^z z 



sm 



(3.10) 



where z = d^z, z' = dpZ and gab{S^) represents the metric over unit 5". We also regulate 
the minimal area surface by stopping it at the UV cut-off z = 5. Now the entropy functional 
becomes 



ICd-3 



4^ 



dpde 



sin^-'^(0) 



yd~l 



{p2 + ^dgzy+p2(d,z)^), (3.11) 



^The cone in d = 4 CFT is also discussed in [23] 
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where is the surface area of the unit {d — 3)-sphere. Note that here we have used 

n = d — 3. From this entropy functional, we can find the equation of motion for z{p, 9). 

= sin(6')z {p^ + i^) z" + sin(6')z (l + z'^) z - 2p^ sm{Q)ziz z 

+id - sin(^) (i^ + (1 + z'^)) + z{d - 3) cos{9)i (i^ + (l + z'^)) 

+pz s\n{e)z {{d - \)z^ + {d- 2)p2 (1 + z'^)) (3.12) 

To proceed further, we use the scahng symmetry of the AdS background. As the background 
geometry is scale invariant and the only scale in our discussion is p, solution of z{p, 9) should 
be of the following form 

z(p,9) = ph{9) , (3.13) 

where h{9) is a function of 9 such that as — )• fi, /i — )• 0. Also at = 0, z achieves its maximum 
value and we have /i(0) = dh/d9\g=Q = 0. We also define /iq such that ho = h{9 = 0). Now 
to extract the cut-off independent term, first we change the integration over 9 to integration 
over h: 

f?-' Js/h,. p A. ftfc"-' ' ^ ' 

Now the plan is to make the integral over h finite. To do so, we find the solution of sin(0) and 
h near the boundary and then see how the integrand diverges in the limit /i — )• 0. Then, we 
subtract the terms up to order 1/h in the integrand of (3.14) to make it finite. At this step, 
we will find that it is convenient to write the asymptotic solution for sin(0) for each d > 4 
separately. Hence, now we work case by case. 

First, we consider d = A and find h and sin(0) in terms of h in the asymptotic limit, where 
/i — )• 0. To do so, we make a change of variable y = s\n.9 and invert (3.12) into equation of 
motion for y = y{h). Using h = y^l — y'^ /y'{h) and h = —{yy'"^ + (1 — y'^)y") /y'^ , we find 
that 

= h{l + h^)y{l- y2) y" - yy' (3 + /j2 + (3 + 5/j2 j_ 2/^4) ^ 2hy^ (l + (l + v'^) 
-h (1 + (1 + h^) y'^) + (3 + h'') y%' - hy^ , (3.15) 

where y' = dy/dh and y" = d'^y/dh?. We can solve this equation of motion perturbatively 
near the boundary and find that 

y = sin(Jl) - \ cosin) cot{n)h^ + f ^(3 - cos(20)) cot^m csc(Jl) log(/i) + yo ) /i^ + Oih^) , 
4 \64 / 

(3.16) 

where we have used y{9 = Vt) = y{h = 0) = sin(r2). In the above expression, yo is a constant 
and its value is fixed by using the fact that y{h) has an extrema sX h = Hq. Further, this 
solution also contains log(/i), which will appear in solutions for even d. As a next step, we 
find h{9) near the boundary in terms of h. For that, we define f{h) = h{9) and write (3.15) 
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as 



0= h{l + h^) yff + /iVl - y^f + (3 + h'') yf 

+h (1 + /i^) - + (3 + 5/1^ + 2/i^) y , 

Now using (3.16) we can solve this equation near the boundary and find 

2tan(Q) 1 



-/i(3 - cos(2f])) csc(2Jl) log(/i) + foh + ..., 



(3.17) 



(3.18) 



h 2 

where /o is a constant. Now using (3.16), (3.18) and d = 4 in integrand of (3.14), we find 
that for small h 



'l + h^ + h^^ + ^°^(^)f»t(^) ^ 



(3.19) 



This implies that the h integral in the entropy functional has only divergences coming from 
first two terms in the limit h — )• 5/ p. We can separate these divergent terms and for d = 4, 
write the entropy functional (3.14) in the following form 

Ici ~ 



5*4 



11 



ih + h 



(3.20) 



where 



H 



dp 



5/p 



dh 



5/ho P Jho 

^ dp f'lP 



sm ( 



1 + /l2 + /l2 + 



sin(r2) cos(J7) cot(r2) 



5/ho P Jho 



dh 



hh3 ' ' ' ' ' ' 
sin(f]) cos(ri) cot(r2) 



8h 



/l3 



8h 



If we series expand Ii in terms of UV cut-off 5, we find that leading term is of order log((5). 
To see that, we use (3.19) and find that in Ii, integration over h is actually finite if we set 
the upper limit h = 0. As all the subleading terms will be of higher order in powers of 5, we 
find that 



log(5) / dh 

ho 



sm ( 



hh^ 



l + h^ + h'^ + 



sin(r2) cos(r2) cot(r2) 



8h 



0{5^). (3.21) 



Simultaneously, we can evaluate I2 and find that 



sin(O) 



4,52 

+ 



1 

16 



cos(Jl)cot(17)log(()7i/)^ 



(3.22) 



^ cos(O) com log(/io) + ] log{S/H) + 0(5°) . 



Now using (3.21) and (3.22) in (3.20), we find the complete structure of divergences in the 
entanglement entropy for cone: 



54: 



4 ^2^3 



ci 



ij2 sm{n) 
452 



1 

16 



cos(O) cot(O) logi5/Hf + 94 logi5/H) + 0(5°) + . .(.3 



•?3) 
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where 



94 



cos(r2) cot(r2) log(/io) 



+ 



dh 



1 + /l2 + /i2 + 



sin(r2) cos(r2) cot(il) 



/i3 



8h 



. (3.24) 



So we find that EE for a cone in d = 4 CFT has a double logarithmic term. We can 
notice from expression of I2 that one of the log terms comes from integration over h and then 
second from integration over p. Here, first integration over h or 9 actually brings us close to 
the cut-off on the smooth part of the entangling surface. Further, when second integration 
over p is performed, we approach to the singularity. This idea is consistent with the fact that 
in EE for even dimensions, we get a logarithmically divergent term according to [16]. 

Now, as a next step, we generalize our discussion to cones in higher dimensions. First, 
we calculate EE for cone in d = 5 CFT. In this case the calculations proceeds similar to the 
d = 4 and the complete expression for the EE is given in the appendix B. However, we find 
that that cut-off independent term is the coefficient of log 5 divergence and it is given by 



5' 



q,{n) logi6/H), 



(3.25) 



where 



55 



4cos2(Jl) ^ sin'^in) 



9hn 



+ 



3 hi 



ho 



dh 



sin^(g) 
hh^ 



1 + /i2 + /i2 + 



sin2(Jl) Acos^{n) 



h^ 



9/i2 



(3.26) 



We further draw this universal term in figure 2. There, loglg,^! is plotted as a function of 
log(sinil) for d = 3 and 5. In the limit — )• 0, we see that loglgsl asymptotes to a straight 
line with slope —1. This implies that for small fi. 



95 oc 



1 



Finally for d = 6, we find that the cut-off independent term is 
iog2| _ 47r3L^ 9cos(Jl)cot(17)(31 - cos(2J7)) 



5, 



IC3 



8192 



logi5/Hy 



(3.27) 



(3.28) 



It is straight forward to see from (3.14) that all the even dimensions will produce a log^S 
divergence. For even dimensions, the number of powers of h in the denominator is odd. 
When separating the divergences, similar to (3.19), it will produce 1/h. This term integrated 
over h and then over p, similar to I2, will produce a log^d divergence. Of course these results 
will persist in field theories in curved backgrounds and for dual gravities with higher derivative 
curvatures. We will discuss these examples in sections 4 and 5. 
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Figure 2: (Color online) We have drawn log \qd\ as a function of log(sinri) for d — 3 and 5. In the 
limit 17 — > 0, it is known that log \q^\ — — log(sinrj) + . . . [19]. For d = 5, we observe similar behavior 
and find that for small il, the linear fit is log {q^l = — 1.01og(sinl7) — 2.1 . 



4. EE for extended singularities 

In this section, we will calculate the entanglement entropy for various creases. Obviously 
the crease can have two types of locus of singularities: flat or curved. By studying these 
examples, we will argue that the contribution from the singularity will be non-zero only if 
the locus is even dimensional and curved. 

In the first subsection, we will study creases with a flat locus of singularity. We will see 
that in this case, we don't find any new logarithmic divergence coming from the singularity. 
The creases we will be considering are k x i?*" with arbitrary m and, Cn x R"^ with {n, m} G 
{{1,1}, {1,2}, {1,3}, {2,1}, {2,2}}. 

In second subsection, we will mainly consider the singular geometries of the form k x 5™ 
and Cn X 5™. In these cases, a generic calculation for arbitrary curvature is tedious. So we 
will work in certain limits where curvature of the locus of singularity is very small and show 
how the leading order terms behave in these cases. In the calculations, we will use following 
approach: If the curvature of the locus is related to 1/Ri, then in the limit i?i — )• oo, the 
boundary geometry will become flat. Now we make Ri finite but keep it very large; and 
calculate the leading order corrections to EE in this limit. In this process, we will find 
that leading correction to EE, in a proper normalization, will be of order 0{1/Ri) but not 
0{1/Ri). Now, this approach has one more advantage. As we are going to argue now, for 
some simple background geometries like R^ x S"^, the coefficient of the logarithmic term at 
order 1 / Rf will be the complete contribution from the kink. We will note in our calculations 
that the only dimensionful quantities in the problem will be the UV cut-off 5 and the IR 
cut-off H. As the coefficient of the logarithmic divergence should be dimensionless and if we 
series expand it in terms of l/i?i, the numerator can either he 5 or H. So the new logarithmic 
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contribution will be of the form 



'new 



= R 




+ 



( 




(4.1) 



In this case, the terms with 5 in the numerator are not UV cut-off invariant and terms with H 
in the numerator are not really the contribution from the kink. That is because H is the scale 
of the bulk part of the entangling surface and a term arising from the singularity should be 
independent of it. So we only need to focus on the leading order correction which is precisely 
the contribution from the singularity. This issue will again be discussed when we rigorously 
calculate EE in section 4.2. 

Let us summarize our strategy for the calculations here: for cases where boundary theory 
is on curved background, we use the Fefferman-Graham expansion to find the dual gravity 
near the boundary. For flat boundary, we don't need to go through this step. Then, we 
calculate the entanglement entropy for the kinks and use the Fefferman-Graham expansion 
to separate the logarithmic divergence from the singularities. We will see that first few terms 
of the Fefferman-Graham expansion will be sufficient to find the leading order corrections to 
the logarithmic divergences in the limit Ri ^ oo. To begin with, we will calculate EE for 
some singular geometries with flat locus in the next subsection. 

4.1 EE for singularity with a flat locus 

In this section, we will calculate entanglement entropy for following two types of geometries^: 
k X and c„ x i?™. We will particularly see EE for cases where n = 1,2 and m = 1, 2, 3. 

4.1.1 Crease k x R"" 

In this section, we will work in the dual Einstein's gravity with bulk metric (2.3) and n = 0. 
The crease k x R^ geometry is given hy 9 £ [— f^, 0], p £ [0, oo] and G [— oo, oo]. We pick 
{p,9,x^) as the induced coordinates over the minimal area surface and the radial coordinate 
is z = z{p, 9). Now the induced metric becomes 




) 



h = 



z- 



(4.2) 



and the EE is given by 





(4.3) 



■p 



'^The EE for wedge x ii^ is also calculated in [23]. 
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where i = dgz, z' = dpZ and we have substituted m = d — 2>. We have also imposed an IR 
cut-off such that G [—H/2^ H /2\. From this functional, we can find the equation of motion 
for z(p, 9) and it turns out to be 

= (p^ + i^) z" + pz {l + z'^) z — 2pzzz z 

+i\{d - l)p + 2zz') + p\{d - l)p + zz') (1 + (4.4) 

Now, in the limit h — )• oo, p is the only scale in the problem. Hence, using scaling symmetry 
of AdS, we can argue that 

z = ph{e) . (4.5) 

Using this in (4.3), we find that 

AttL'^-^H'^^^ dp f^''' „ 1 



where we have put the UV cut-off at z = 5 and defined /iq = h{0). Further, we have used 
(4.5) to find the limits of the integrations and the extra factor of two comes from changing 
the range of the integration from 9 G [— f^, 0,] to 9 & [0, Q]. Using (4.5) in (4.4), we find the 
equation of motion for h 

h{l + h^)h + {d- l)h'^ + + dh^ + (d-l) = 0. (4.7) 

In this case, although it is not explicit, we can find a constant along the 9 translation from 
(4.7). It is straight forward to see that 

is a conserved quantity and it satisfies the equation of motion. Now we can use (4.8) to 
separate the divergences in the EE (4.6). Using the fact that h decreases near the boundary 
and hence h should be negative, we find 



^/rT7?./(l + /l2)rf-2 _ /l2(rf-l) 



Using this in integrand of (4.6), in the limit /i — )• 0, we find that 



(4,9) 



uud-^ hd-l 9 " ^ ' 



Now in (4.6), we can add and subtract l/h'^~^ in the integrand and write 



(4.10) 



where 



h 



H 



dp 



S/ho P'^ ^ J ho 



5/p 



dh 



^/h? + /i2 + 1 ^ 1 



In (4.11), first few term are coming from h and p integration of l/h'^^^ . To find the divergences 
in /i, first we represent the integrand by J{h) and from (4.10), J{h) ~ h'^~^ as /i — t- 0. We 
perform the integration by parts and write Ii as 



{d - 3)Hd~^ 
1 

{d - 3)Hd~^ 



S/H ^ 

dh J{h) — — — — 
ho (" - 3) J s/ho P' 

dh J{h) - — — ^ h 



ho 



{d-3) 



(4.12) 



Now to separate the divergences in I2, we make a change of variable from p to q = 5/p and 
then Taylor expand the terms around 5 = 0: 



1 f/w 

^ / dqq'^-^Jiq) 
J ho 

^ j\qq''-^J{q) + ^{q''-^J{q) 

^ f dqq''-^J{q) + {0){5^). 
J ho 



q=5/H 



(4.13) 



In the above expression, the integral over q is finite because from (4.10), q'^ ^J{q) ~ q"^"^ ^ 
for small q. Combining (4.11)-(4.13), we can write 



d—\ fjd—3 



47r L^^-'H 



H 



{d - 2)5^-^ 



+ 



.,d~3 



ho 



1 

ho 



{d - 3)5'^-3 



+ 0(<5°) 



(4.14) 

Remarkably, we find that as soon as we add a flat locus to the kink, the log divergence 
disappears. However, there is a new divergent term of order l/5'^~^ which does not arise in 
smooth entangling surfaces. This clearly shows that logarithmic contribution from the crease 
depends on the curvature of the locus of the singularity. Further, there can be a logarithmic 
term from the smooth part of the entangling surface in even d. However, the entangling 
surface in A; x is flat everywhere apart from the singularity, and hence such contributions 
vanish. This example clearly shows that when the locus is flat, there is no contribution 
from the singularity. However, we will see ahead that as soon as we turn on the curvature 
of the locus, the logarithmic divergence will appear only for the cases where locus is even 
dimensional. 

4.1.2 Crease c„ x R"" 

In this section, we will calculate EE for the geometries c„ x i?™. Although such a calculation 
can be extended to arbitrary n and m, we will particularly focus on n = 1, 2 and m = 1, 2. To 
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begin with, we consider the metric (2.3) for the dual gravity with arbitrary m. The conically 
singular geometry Cn x i?™ is given hj p £ [0, oo] and 9 G [0, f]] in gravitational dual (2.3). We 
assume that the induced coordinates over the minimal area surface are (p, 9, x^), where .^*'s 
are coordinates over unit sphere and x-^'s are on ii™. As we have rotational symmetry over 
the sphere S"", the radial coordinate for the minimal area surface will be given hy z = z{p, 9). 
Now, the induced metric over the surface is 



(1 + .-) 

L2 



-z'i 



L2 

-71 



72 



L2 



(4.15) 



where gabiS"^) represents the metric over the unit sphere S'" and the EE is given by 

2-K L^^-^H "^nr 



if"' 



dpd9 ^"7_7^V (^^+p^(i+^^^)) 



(4.16) 



where Vtn is the area of the unit n-sphere, z = dgz and z' = dpZ. Note that we have integrated 
over the x*'s and used the IR cut-off G [—H/2,H/2]. Now for this case, the equation of 
motion becomes 

= (p2 + z" + p^z (1 + z - 2p^ziz'i' + {d- l)p^ (i^ + p2 (1 + z'"")) (4.17) 
+z {d' cot{9)i^ + {d' + 2)pz^z' + d! p^ cot(0)i (l + z'^) + (d' + \)p^z' (l + z'^)) . 



Once again, here we can use the scaling arguments and say that 

z{p,9) = ph{9). 

Then, EE reduces to 



2tt L'^^^m'^n 



H 



dp 



lc„x_R'" 

and equation of motion for h becomes 



S/hti P J ho 



^Ip „ sin"(^)vA2 + /i2 + i 

dh — 



1 



(4.18) 



(4.19) 



= h{l + h^)h + ncot{9)hh^ + {d + nh^ -l)h^ + ncoi{9)h{l + h^)h 

+{n + l)/i^ + {d + n)h^ + d-l. 



(4.20) 



In (4.19), we have changed the integration from 9 to over h. We have also introduced the UV 
cut-off z = 5 and defined /iq = /i(0). 
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Now we set n = 1 and d = 5 (that also mean that m = 1), and calculate EE for the 
singular surface ci x R^. First, we need to find y = sm(6) near the boundary in terms of h. 
For that, we invert the equation of motion (4.17) and get 

= h{l + h^)y{l- y2) y" _ 2 (2 + 3/i2 + h^) yy'^ (4.21) 

-h (1 + h^) (1 - 2y') _ (4 + h') y (l - y^) y' -h{-l + y'f . 
Now solving this equation perturbatively near the boundary, we get the solution 

y = sin(f]) -hi'^ cos(O) cot(J7) - |^^''(19 - 5 cos(20)) cot^iVt) csc(O) + 0{h^) , (4.22) 

where we have used that at /i = 0, y(0) = sin(r2). Further, we can define h{9) = f{h) and 
write the equation (4.20) in the form 

Q = h{l + h^) yff + hy'l-y^f + (4 + h^) yf 

+h (1 + y/l -y2f + 2[2 + 3h^ + h^)y, (4.23) 

where y' = dy/dh. Using y from (4.22) and solving this equation near the asymptotic bound- 
ary, we find 

/ = -^^^^^ + J/i(8-cos(2f]))csc(2J7) + /o/i2 

3_/i3('435 _ cos(217) + 52 cos(4Jl)) csc^(f^) sec(17) + 0(/i^) , (4.24) 

216 

where /o is a constant which is fixed by the condition /(/lo) = 0. Now it is straight forward 
to see that near the boundary 



sin(^)Vl + /i2 + /i2 sin(O) cos(J7)cot(17) 1 
and we can use it to make h integral in entropy functional finite. So we write EE as 



2 , 



H^sm{n) cos(^^) cot(Jl) , 




2/igcos(Q)cot(J7) (1 - log(/io)) - 9sin(17) 
^ 18/^p 

sm^'{eWh^ + h^ + l sin(O) cos(0)cot(0) \ ^1 

where we have already performed the integrations over some terms. Now let us call the term 
with integration Ii and the integrand J^{h). Then, near the boundary J5(/i) ~ 0(71°). Now 
using integration by parts, we can write 

h = dhMh) 

Js/ho P Jho 

/ dhMh)-5 -iMh)\h=5/p 
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We further make the coordinate transformation p = and Taylor expand the second term 
in 5: 



h = -T? / dhJ^{h) + - I dqqh{q) 
-H Jho (> Jho 

= 1^ dhMh) + ^ l^qq^Mq)- j^cos{n)cot{n) + 0{6), (4.27) 

where we have used (4.25) to get the final expression. Combining (4.26) and (4.27), we can 
write 



55' 



H'^s\n{Q) 1 /"^ cos(l^)cot(0) log((5/iJ) 

dqqJ^lq) ' 



2hl cos{Vl) coi{n) (1 - log(/io)) - 9 sin(^7) . 



(4.28) 



Note that in the above expression, we have a new divergence of the form \og[5/H)/5 which 
does not arise in EE for smooth entangling surfaces. This term should be a contribution 
from the singularity. Further, we note that as soon as we add a one-dimensional locus to the 
conical singularity, both double log and logarithmic divergences disappear. Recall that in the 
previous case for cone ci in d = 4, i.e., eqn. (3.23), we got a double log and log terms. 

We can easily generalize the above calculations to the crease ci x E?. For this case, the 
integrand near the boundary and complete expression for EE are given by (B.5) and (B.6). 
Once again we find a divergent term of the order \og{5/H)/5 in EE. Also, now the cut-off 
independent term is 

log, 47r2 L^H^ 3(13 - 19 cos(20)) cot^(^^) csc(O) ^ , 
^6 l^.R^ = ^5 ^I^^^^ ^og{5/H). (4.29) 

Note that for singular geometry ci x i?^ in d = 6, we only have log divergence instead of 
log^J as compared to the cone ci in equations (3.23). Also, recall that in (3.28), we saw 
that a cone without a locus, C3 in d = 6 gives a log^(5 divergence. In (4.29), this log^(5 term 
disappears because of flat locus we have added. We will see in section 4.2.2 that as soon as 
the curvature of the locus turned on, the log^(5 divergence will reappear in (4.147) in d = 6. 
In (4.29), the logarithmic contribution can be attributed to the fact that in even dimensions, 
the entangling surface has a logarithmic divergence. We can explicitly verify that logarithmic 
term in (4.29) comes from integration over 1/h in expansion (B.5). As p need not to be near 
the singularity, this contribution is from smooth part of the entangling surface. Further, it 
is straight forward to find that for cone ci x B? \n d = 1 dimensional CFT, the EE has no 
double log or log divergence. Which is expected because for entangling surface in odd d, there 
is no logarithmic divergence and as the locus is odd dimensional and flat, there should not 
be any logarithmic contribution from the singularity either. 
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We can further extend these calculations for the case n = 2 and m = 1,2. For crease 
C2 X R^, we find that the universal term in EE is a logarithmic divergence and it is given by 

log, 8 TT^ L^H (7 - 9 cos(20)) cot\n) 
^6 = ^5 ^ ^og{6/H). (4.30) 

We also find a new contribution from the singularity at the order 1/6. In x R^, the locus of 
the singularity is odd dimensional and hence singularity doesn't contribute in the universal 
term. However, smooth part of the entangling surface contributes through a log in d = 6 and 
this is the contribution (4.30). It can be verified that this log term arise from integration 
over h away from the singularity. Further, we can also calculate EE for crease C2 x R^ and 
find that there is a new divergence from singularity of order 1/5^. However, there is neither 
a log nor a log^ term in EE. Note that for crease C2X R?, the locus of the singularity is even 
dimensional but fiat. So singularity does not contribute through a log term. However, as we 
will see in (4.152) that for a curved locus, we get a log contribution. 

4.2 EE for singularity with curved locus 

In this subsection, we will consider several singular embeddings which have a curved locus. As 
soon as the curvature of the locus is turned on, the double log and log terms from singularity 
will make appearance. In this section, we will consider the creases A; x S and c„ x S, where 
locus S will take the form or 5*""^ x RP. These cases will be slightly more involved 
and hence we will do two things: first, we will always work in the limit where curvature of 
the locus is very small and we will do the calculations perturbatively. In certain cases, we 
will see that these perturbative calculations are sufficient to pick the complete contributions 
from the singularity. Second, for calculations in this section, we will foliate our minimal area 
surface in different way. In all the previous cases, the induced coordinates on the minimal 
area surface were {p,9, . . .) and we assumed that the radial coordinate z = z{p, 6). However, 
we will find that it is more convenient to do the calculations in a coordinate system where 
bulk radial coordinate z is one of the induced coordinates on the minimal area surface and we 
have p = p{z, 9). We have shown in figure 3, how these different induced coordinates foliate 
the minimal area surface in the bulk. Note that the new set of induced coordinate (z, 0, . . . ) 
are not the well-defined coordinates on the boundary as p = />(z, 9) will be multivalued for 
z = and 9 = ibf]. However, until unless we put the UV cut-off, we find that p = p{z, 9) is 
a well-behaved function and we can work with it. Now, in the next subsection we will begin 
with the geometries A; x S and then we will move on to conical singularities c„, x S. 

4.2.1 Crease A; x S 

In this section, we will mainly consider the geometries k x S'^, k x R x S"^ and k x S^. We 
will see that singularities with even dimensional locus will contribute through a logarithmic 
term. 

To begin with, let us consider d = 5 CFT on background R'^ x S^. Before we construct 
the singular entangling surface in this geometry and calculate holographic EE, we need to 
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Figure 3: Panel (a) shows how minimal area surface is foliated when we have the induced coordinates 
{p,9, . . .) and minimal area surface is given by z = z{p, 9). Panel (b) shows the foliation of the minimal 
area surface when induced coordinates are (z,9, . . .) and we have p = p{z, 9). Note that z = z.^ is the 
maximum value of z on the surface such that p{0, Zm) — H. 



find the dual gravity. So we begin with the action for six-dimensional dual Einstein gravity 



1 



(4.31) 



for which, the equation of motion is given by 



Rt^v - ^gt^u = . (4.32) 

Now we ansatz that the bulk metric, which has the boundary x S^, is of the form 

ds"^ = — {dz'^ + /i(z) {dt^ + + p'^dO'^) + f2{z)Rl d^l) , (4.33) 

where d^\ = d^Q + s\v?{^Q)d^\ represents the metric over two-sphere and, /i and /2 are 
functions of the radial coordinate. Here Ri is the radius of the sphere on the boundary and 
in the limit i?i — )■ c«, we recover the flat boundary. To find /i and /2, we use the Fefferman- 
Graham expansion near the boundary z = 0. The idea is to insert the ansatz (4.33) in the 
equation of motion (4.32) and then solve it near z = to find 



z^ 17z^ 
12Rl ^ 576Rt 

z2 



f 1 = 1 H wT H + . . . and 



Knowing the dual geometry near the boundary, we can calculate EE for the geometry kx S'^. 
On the boundary, this geometry is defined by p G [0,H] and 9 £ where H is the 
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IR cut-off. Now to define tlie minimal area surface, we use sliglitly different coordinates 
compared to the previous sections. We choose the induced coordinates to be iz,6, ^o, ^i) and 
because of the rotational symmetry on the sphere, p = p{z, 9). Although, p{z, 9) is not a well 
defined function at the boundary z = but we can definitely work with it until we impose 
the UV cut-off and assume that the boundary is at some finite z = S. With this coordinate 
choice, the induced metric becomes 



h 



ihp" + 1) 



' ip' + P') 



,2 



■sin^(^o) 

where p = dgp^ p' = dzP and hence, the entanglement entropy is given by 

/iT/2 



(4.35) 



dz d9 



(4.36) 



We can easily find the equation of motion of p{z, 9) to be 



= 2zhhp {p" + (?) p" + 2zhp (1 + /ip'2) p - Azhhppp'p' + 2zhfW {f? + p2 (1 + f^p'^)) (4.37) 
+/2 {-Azf? + p (-8/1 + 3z/0 f?p' - 2zp^ (1 + + p^p' (3^/1 _ 8/2p'2 + (-4 + zf[p'^))) , 

where p = dgp, p" = dip and p' = dzdgp. 

As discussed earlier, we will be working in the approximation Ri — t- oo. In this approx- 
imation, we will find that leading order correction to the EE will be of the order 0(l/i?f). 
To show this, first we need to find the form of the solution p{z, 9) in this approximation. For 
that, we make following ansatz: 



(4.38) 



This ansatz is made on the basis of the available dimensionful quantities in the problem. Now 
we are going to argue that in the following expansion, the 9 dependent functions g2n are zero. 
To see that, we insert the ansatz (4.38) in the equation of motion for p{z, 9) and find following 
equations of motion for /i, g2 and g^,: 



h{l + h^)h + Al? + (1 + h^) (4 + h^) 
{l + h^f g2+2h{7 + 2h^) {l + h^) hg'^ 

+g2 (lAh^ - 2/1^(8 - h^) -h^ - 31h^ 
(1 + h^)^g3 +Ah{A + h^) (1 + h^) hgs 



0. 



16' 



-93 ( /i' + /i' ( 35 



2/i^ 4 5 - 6/i^ 2h 



Si 



(4.39) 
(4.40) 

(4.41) 
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where Si is given by 



+ (^32 - 15^2) +h^(8 + + 2/1^ (15 + + 2 (5 + SA^)) . (4.42) 

Note that here we have arranged the equations of motion such that on the left hand side, we 
have the homogeneous part of the equation and on the right hand side, we have the source 
terms. So Si is source for and in (4.42), terms in first hue come from the corrections to 
fi and /2, and terms in last two lines are second order in g2- We also note that in (4.40) 
and (4.41), we have used equations of motion to eliminate the second derivatives of h and 
g2. Now we notice following points: first, the corrections to fi and /2 are of even powers of 
and hence they only source g2n+is. Second, g2n are sourced only by 521 's, where i < n. 
Third, in flat boundary, precisely in the limit Ri — )■ 00, we have h ^ and 52 = 0. So above 
arguments conclude that as we make Ri finite, only (72n+i(^)'s will be turned on and all the 
g2n = 0, which clearly are solution of these equations of motion. 

Now to separate the logarithmic divergence in entanglement entropy, first let us write 
p = po{z, 9) + pi{z, 9)/Ri, where po = z/h{9) and pi is higher order corrections for large Ri. 
Using this and (4.34) in (4.36), and keeping only the leading order terms in we find that 

(po {5z^po - 24pi) - 24popi (l + p'^^) + pj (5^2 + iz^p'^'' - WoP'v\)l . 

2Az^RWpI + pU1 + {p'o)') ™ 



m 



where 5 is the UV cut-off and e = e{z) is defined such that p{z, 0, — e) = H and further, z. 
is defined such that p{zm, 0) = H or e{zm) = We can also insert this ansatz in (4.37) and 
find the equation of motion for po and pi by series expanding it in terms of Ri . In the second 
term of (4.43), we can convert p'l and pi into pi using the integration by parts and then use 
equation of motion for po to simplify the coefficient of pi to zero. This process will leave us 
with some boundary terms and we find that 



^^IkxS^ = M / / 

P J Zm J —fl+e 



Vp^+p2(l + (p'o)2) (5p2 + p2 (5 + 4^/2) ) 



Z 



4 



24i?2zVpg + P^(l + (/^'o)') 



I Popi \ _^ d ( PoP'oPi 



09 \i?2^4^^2 + ^2(i + (^.)2) / dz U2^Vpg + pg(l + p'^2) 



(4.44) 
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Now we can insert the ansatz pQ = z/h{6) and pi = z'^g3{6) and simplify the functional to 



^^"^ J Zm J ho h 



^^^^ 2ARlzh'^yJl + h? + h? 



(4.45) 



where we have changed the integration limits from (— il, fi) to (0, i7) and also changed the 
integration variable to h{0). Note that the boundary term with derivative with respect to z 
turns out to zero up to leading order. We have also defined /iq = /i(0) and /iic(p) = h{VL — e) 
and used /i'i(O) = in getting the boundary terms. 

Now to separate the logarithmic divergences, we consider all the contributing factors 
one by one. We will see that the first term in the first line of (4.45) will not contain any 
logarithmic divergence. The second term and the boundary terms will contain logarithmic 
divergence. Note that we also need to be careful and consider the divergences coming from 
the limits of the integrals. To begin with, we first study the behavior of h and g^, near the 
asymptotic boundary. The equations of motion for h and 53 are given by (4.39) and (4.41) 
with g2 = 0. Now similar to (4.8), we find that equation of motion for h can be integrated 
once to get 

K,.^±£L, (4.46) 

where is a constant and it can be further related to /i(0) using /i(0) = 0. 

To extract the logarithmic divergence, we will only need the asymptotic behavior of h 
and 33. Hence we solve g'3 in terms of h in asymptotic limit, where h is small. To do so, we use 
(4.41) and change the variable from to /i by expressing 53(0) = d'^g^/dO'^ and 53 (^) = dg^/dO 
in terms of g3{h) = d^g^/dh? and gz{hl) = dg^/dh. We further use (4.46) to express hl'{6) 
in terms of and h. Finally, the equation of motion for g3{h) becomes 

= 12h^ (1 + /i^) 2 [h^Kl - (1 + h^) 3) 53 + 12h^ {h^ (1 + h^) (16 + 5h^) Kl - 4 (l + /i^) ^ (3 + h^)) 
-12h (2 (1 + h^) 3 (12 + h^ + h^) - (44 + 17h^ + 3h^) 53 

+ (1 + h^) 3 (3 + 7h^) + hh^Kl . (4.47) 

Now this equation has two solutions when we solve it perturbatively in the limit /i — t- 0. The 
leading terms of these solutions go like and However, 53 must be such that p is 

finite in the limit ^ — ?• and (5 — ?■ 0. As g^ appears at order 5^ in (4.34), it can only go like 
Hence, one of the constants, which is the coefficient of l//i®, is fixed to zero. As a 
result, the asymptotic solution turns out to be 

63 1 + 8863 ^ 4 + 7263 4 + 7263 ,3 , UA9l\ 

^^ = ^ + ^6)^ + ^89-^-^93-^ +•••' ^'-^^^ 
where 63 is a constant such that gz{0) has an extrema at = 0. 
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Before we begin discussing divergences of various terms, we find the series expansion of 
/lie in terms of the UV cut-off 6. As hic = h{Q — e) and z = 6 at the cut-off, we can use 
(4.38) and perturbative solution (4.48) to find following series expansion for hic 



hicid) 



1 , b3H\ 
Rl) 



6 + 



(1 + 



56HRI 



(4.49) 



Note that we have kept only leading corrections in Ri at any order in 5. 

Now we return to (4.45) and analyze the divergences for each term. First, we use (4.46) 
in the integrand of first two terms of (4.45) and find that in the asymptotic limit 



5/i2 + 5/j2 + 4 



1 



1 



2ihh^Vl + W+¥ 24/i2 16 

So we can make the integrands finite by organizing the terms in following form 



(4.50) 
(4.51) 



h 



dz 
dz 



hi 



dh- 



ho 
ho 



hh^ 



dh 



Vl + h^ + 1 



+ 



dz 



1 

hic 



1 

ho 



A + 1'2 



2 ) 



(4.52) 
(4.53) 



and 



Ik + I'A 



dz 


rhic 


z 


J ho 


dz 


/•hic 




'ho 



dh 



dh 



bh^ 5/i2 + 4 

2Ahh?^Jl + h? + h? 
5A2 + 5/i2 + 4 
24A/1V1 + + /i2 ' 24/i2 



+ 



+ 



dz 



1 

2Az \h^c 



1 

ho 



(4.54) 
(4.55) 



In (4.53), /( and I2 represent the first and second integrals in (4.52). Similarly in (4.55), 13 
and I'^ are the first and second integrals in (4.54). Now first we consider /[. We differentiate it 
with respect to the UV cut-off 5 and look for 1/5 divergent terms. After taking the derivative, 
we find 



dl[ 

l6 



1 

P 
1 

P 
1 



hiciS) 



ho 


I 

ho 


ho 



Vl + h^ + h? 1 
A/l2 + 



dh 



dh 



dh 

'^l + h? + h? 



+ 



+ 



1 
1 

h? 



+ 



1 dhU5) 
52 d6 



+ 



Vl + /l2 + h? 1 



- h=hic{5) 



2m \ m ' Ri 



5* + C'((5*' 



(4.56) 



-27- 



where we have Taylor expanded the integrand in the second hne and in the third hne, we 
have used (4.46) and (4.49) to find the leading order divergence. Here we don't have any 
term of order 1/5 and of \og{5)/5, which come consecutively from divergences of order \og{5) 
and log((5)2 in I[. 

Similar to /(, we can take a derivative of with respect to 5 and use (4.49) to find 



1 



1 



1 + 8863 H 4 + 7263 



(4.57) 



d5 y Rl J 54 ho5^ 56Rj 189 RjH 

So we find that Ii doesn't have any logarithmic divergence. We can use similar steps to find 



dl'3 


_ 1 


d5 


^ J ho 


dli 


5H 


d6 


2452 



dh 



1 



5A2 + 5/i2 + 4 5 
+ 



24/l/l2Vl + /l2 + /i2 24/i2 



^ + ] 6^ + 0(<5^|4.58) 



i/5 \16H^ IGRj 



Rl J 24/10 <5 

We can also use the same procedure on the boundary term and find that 



(4.59) 



d 
d6 



dz 



Vl + /l2 + h'- 



bEl (^ _ 3^3^' \ (1 + 8863) H {-4b3H^ + Rj) 
54 Rj 56Rj6^ 

4 + 7263 (27 + I663 (521 + I7IOO63)) H 



189H 84672i?2 
So from (4.56) - (4.60), we can find the logarithmic divergence in the EE for kx S'^ geometry: 



5; 



log\ 



16 7T^L^ 



4 



dh 



ho 



5^2 + 5/i2 + 4 



+ 



24^2^1 + /l2 + /i2 24/i2 



+ 



24/10 



log(5) . (4.61) 



Prom (4.56) - (4.60), we also notice a divergence of order 1/6'^ in EE which does not appear 
in EE for smooth entangling surfaces. Note that such a term also appeared in EE for k x i?"* 
in (4.14). We further notice that logarithmic contribution in (4.61) is of next to the leading 
order for large Ri. As we are working in odd dimensional spacetime, there is no logarithmic 
contribution from the surface itself. 

Now we turn to our next example. Having seen the logarithmic contribution from a 
even dimensional curved locus, now we consider odd dimensional locus. We will calculate the 
entanglement entropy for k x geometry in CFT on R^ x S^. For this case, the locus is 
and we will see that there will be no log contribution from the singularity. However, as 
the CFT is in even dimensional spacetime, we should be getting a logarithmic contribution 
coming from the entangling surface. In this case, the metric for the dual geometry is given 
by 

ds^ = {dz^ + fi{z) {dt^ + d/>2 + p^d9^) + f2{z)Rl dnfj , (4.62) 
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where dfig = d^^ + sin^(^o)'^^i + sin^(Co) sin^(Ci)f^C2 the unit three-sphere and /i and /2 
are following 



/i 

/2 



69z^ 

1 + TTTTT^ + T^T^TT^ + • ■ ■ , 



1 



20Rl 1600i?f 
20Rl ~ 1600i?f 



+ ... . 



(4.63) 



Once again we choose the induced coordinates on the minimal area surface to be {z, 6, (,o,S,i,S,2) 
and assume p = p{z, 6). Then the induced metric over the minimal area surface will be given 
by 



h 



S ihp'^ + 1) ^PP' 



'-^PP' 



"^{P' + P') 



(4.64) 



where p = dep, p' = dzP and now the entanglement entropy is given by 



5, 



6|fcx53 



1 



(4.65) 



(4.66) 



We can easily find the equation of motion for p{z, 6) to be 
= 2zhhp (p2 + ^2) p" + 2/2P {z + z/ip'2) p - Azhf2ppp'p' 

+3zfipf^p' (p2 + p2 (1 + _ f^(^4,p2 + ^ (lo/i - 3z/0 p2p' + + ^_^^^/2) 

-/^^p' {3zf[ - lO/iV'2 + 2/i (-5 + zf[p'^)) 



where p = dgp, p" = dip and p' = dzdep. Once again, we can write p = po + Pi/Ri and keep 
only the leading order correction to the entanglement entropy. We can further series expand 
the equation of motion (4.66) to get equation of motion for po and pi. We find that, using 
these equations of motion and integration by parts, the entropy functional can be simplified 
to 



At^'-^L^RI 



XS3 



£5 



+ 



dz 



de 



d 



PoPi 



V/3^ + p2(l + p^2) 3 (6p2 +pI{q + 5p^2)) 



40i?2,3^^2 + p2 (1 + ^/2) 



+ 



d 



pIp'gPi 



de yRiz^>^p2 + p2(i + p'^2)j dz yRiz^^pi2 + p2{i + p'^2) 



^4.67) 



Further, we can substitute po = z/h{e) and pi = z^g^ie). The equations of motion for h and 
5(3 from (4.66) becomes 

h{l + h^)h + bh^ + (1 + h^) (5 + /i^) = , (4.68) 
(1 + /i^) 2(^3 + 2/i (9 + ll/i2 + 2/1"^) 

-53 (25 + 45/1^ + 21/1^ + /i'^ - (27 - /i^ + 2/i^) lA = Si , (4.69) 
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where source terms for 53 are 

3 f 10 + 19/i2 + 9h^ + (4 + 9/i2) hA 
Si = , ^ . (4.70) 

Here, as in (4.8), we can integrate the equation of motion for h once and write it as 

= \ ' . , (4.71) 
W'^/l + h'^ + h? 

where is a constant and we can relate it to /iq by using /i(0) = 0. We can further simphfy 
(4.67) to 



where 



f /■"■• dh v'l + feZ + ft2 



h = - dz ^ ^ ^ . , (4.74) 

Ao /l 40i??z2/l2Vl + /l2 + ^2 



3 



1' 

dz gsh 



(4.75) 



n- 



/ ho 



Note that the last term I4 is not being integrated over z. That's because hic{z) appears in 
the limit of 6 integral. Now, we can use (4.71) in (4.69) to find following equation of motion 
for (73 in terms of h: 

= 20/1^ (1 + h^) 2 {h^^Kj - (1 + h^) ^) gs + 20h\l + h^) (h^° (I8 + 5h^) Kj 

- (1 + h^) ^ (13 + 4/i2) )^3 + 20/i (/i^° (52 + l^h^ + 3/i^) - (l + h^) ^ (27 - /i^ + 2/i^)) 93 
+3 (1 + 4 (4 + 9/i2) + is/jio^l _ ^4 

This equation can be solved perturbatively for small h near the asymptotic boundary. We 
find that this second order equation of motion has two different solutions which go like 
and 1/h?. As /> is finite at the boundary, the solution can go only with power 1/h?. So the 
solution near the boundary becomes 

^63 3 + 140^3 ^ 3 + 3663, _ 1 + 1263 ,3 + 5 + 6O63,, 
' ¥ 80h 64 96 1536 ^ ^ 
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where constant 63 is fixed by the condition that gsiho) = 0. We use this in p = z/h + z^g^/ 
and evaluate \i at z = 5, p = H and h = hic{6). By inverting the relation and keeping only 
leading order terms in Ri, we find 

Now we use this to separate logarithmic divergence. Using (4.71), we find that the integrands 
in Ii, I2, I3 and I4 are of following form for small h: 

Vl + h^ + h^ 1 1,-2,8 , .o..9^ 



hh^ 2 
3 (5 + 6/i2 + QhA g 3 



40A/iVl + + A2 20/i2 20 

'-^ ~ + + A (1 + 12^3) /,_!(! + 1263) /,3 + o(/,5) 

/i'^ 80/i 64 ^ ^ 96 ^ ' ^ \ I 



93 



, = ~ bsKlh^ + 0{h'^) . (4.80) 

hhVl + h^ + h^ 

Using these relations, now we can make the integral over h finite and organize the terms in 
following form 




1 + /i2 + A2 _^ 1 \ /"^ / 1 1 



h'^J J^^ z^ \hiciz) ho 

(4.81) 

40/i2\/l + /i2 + /i2 ' 20/i2j 20 J,^ z^ \hiciz) ho^ 

(4.82) 

In (4.81), I[ and I2 represent the first and second integral in the expression of Ii. Similarly 
in (4.82), /g and I'^ are the first and second integral in expression for 12- We can further take 
derivative of /"s with respect to 6 and then Taylor expand the integrals to get 

dl[ _ 1 f^yi + h-^ + h^ Kl[Qh;H^ + Rl) 



dik - b3H^\ 1 (3 + 14063)g 3 + 36fc3 , .^.,0. 

dS Rl ) hod^ 80 Rl5^ UHRl5^ ^ ^ ' 



d6 6^ Jh, \^40/i2Vl + /i2 + A2 j 20Hml 

dl!i / hH^\ 9 9 (3 + 140 63) 



dS - 20 5-^ V RlJ^^^^ imORl5 + • 
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Also, we can simplify the boundary terms in a similar fashion and find that 
~ + 1 ' 03 + 



( 3 + 36 63 (9 + I44O63 + 26800 61) g \ 1 
+ 1^~64^^ 6400i?^ J ^ + ) ' ^^-^^^ 

§ = -1 r dh ^ /\ . + 0(<^^) . (4.88) 

Now using (4.72) and (4.81)-(4.88), first we notice that there are new divergent terms of order 
1/5^ and 1/5 which does not arise in EE for smooth entangling surfaces in d = Q. We also 
find that the logarithmic divergence of the entanglement entropy are 

fr"L..3 = '^f^(o + o(ij))log(<S). (4.89) 

Remarkably, the contribution to the logarithmic divergence at the leading order is zero. At 
higher order, i.e., at order O (l/i?f), there should be logarithmic contribution coming from 
the bulk part of the entangling surface. According to [22], the logarithmic contribution from 
four-dimensional entangling surface will have the coefficient of the form J dx'^C, where C is 
a combination of various curvature squared terms. These terms will contribute with l/Rf 
and so does the smooth part of the entangling surface. So up to leading order, there is 
no logarithmic contribution from the singularity because the locus is odd dimensional. We 
might ask if there can be any contribution from the singularity alone at higher order in 1/Ri. 
However, we discard such a possibility. It is easy to see that all the higher order terms will 
be of even powers of 1/Ri and to make the coefficient dimensionless, the only other available 
scale will be H. As H is related to the size of the surface, all such contributions actually 
result from the contribution from the bulk surface and not the singularity alone. Hence, we 
expect that there is no contribution from the singularity in this case. 

Next, we see one more example where odd dimensional locus doesn't contribute to log- 
arithmic dimensions, in spite of non-zero curvature. We are going to consider k x R^ x 
geometry. We consider the background geometry for CFT to be R^ x S'^. Then the bulk 
metric for boundary i?^ x S"^ is given by 

ds^ = — {dz^ + /i(z) {dt^ + dp^ + p'^de'^ + dx"^) + f2{z)Rl dfl^) , (4.90) 
where d^l2 is line element over 5^ and now /i, /2 become 

^ 20Rl ^ lOORj ^ ■ ■ ■ ' 

z"^ Iz^ 

/2 =1 ^ - + . . . . (4.91) 

^ 5Rj 800i?^ ^ ^ 
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The k X X S'^ geometry is given by G [—J], J]], x £ [—00,00] and p £ [0,oo]. We put 
IR cut-offs on X and p directions sucli that x € [—H/2,H/2] and p £ [pm,H], where pm is 
related to the UV cut-off 5. We choose (z, 0, x, Ci) as induced coordinates on the minimal 
area surface with p = p{z,9). Then, the induced metric over the minimal area surface will 
become 



[P + P ) 



where p = dgp, p' = dzP and now the entanglement entropy is given by 



(4.92) 



kxR^xS^ 



4 



I dzd9^Vp^ + PH^ + fiP'^)- (4.93) 



Note that here we have already performed integration over x. In this case, the equation of 
motion for p{z, 9) becomes 

= 2z/i/2p (p2 + p2) p" + 2zf2p (1 + fip'^) p - 2zf2 + 2p2) p' 



(4.94) 



+2p {2zM[ + h (-5/2 + zf^)) (p2 + p2) pi _ 2z/i/2pV" + /ip' (3^/2/( + 2/1 (-5/2 + zf^)) p 

where p = d^p, p" = dip and p' = dzdep. Once again, we can write p = po + pi/Rf and keep 
only the leading order correction to the entanglement entropy. We can further series expand 
the equation of motion (4.94) to get equation of motion for pQ and pi. We find that, using 
these equations of motion and integration by parts, the entropy functional can be simplified 
to 



/3 



Sa 



Stt'^L^RIH 



dz 



+ 



PQPl 



d9 



/p^ + p^(l + p'o2) 3 (6p2 + pI{Q + Sp'o^)) 



+ 



d 



mRlz^^pi + piii + p',-^) 

(4.95) 



pIp'gPi 



39 \R2,^^p2 + p2^l + p'^2^ dz \R2,^^p2+p2(l+p'^2 



Further, we can substitute po = z/h{6) and pi = z^g^(9). The equations of motion for h and 
5(3 from (4.94) becomes 



h{l + h^)h + hh^ + (1 + h^) (5 + /i^) = , 



(1 + h^) ^gs + 2h{9 + llh^ + 2h'^) hg. 



-g-i (2b + 45/1^ + 2lh^ + - (27 - + 2/i^) 
where source terms for 53 are 



Si 



(10 + 19/i2 + 9/1^ + (4 + 9/i2) (^h 



20h 



(4.96) 
(4.97) 

(4.98) 
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Note that homogeneous part of the equations of motion for h and are same as the case of 
kink in x 5'^. Similar to (4.8), we can further integrate the equation of motion for h to 
get a conserved quantity 

Kg = \ ' . . (4.99) 
The simphfied expression for the entanglement entropy then becomes 



where 



/.= (4.01) 



ho h z^K^ 



s /./lie ( 5 + 6/1^ + 6/i2 
dz ^ ^ / - . (4.102) 



/3 



dz g-^h 



dz f^^"" dh 53 



(4.103) 
(4.104) 



J ho h h^/l + h'^ + h? ' 

Now we can convert the equation of motion for g^{9) into equation of motion for 53 (/i) and 
solve it perturbatively near the boundary in the limit /i — )■ 0. We find that 

63 I + I4O63 I , I , , , o {h + imhz)h^ , 

03 = -4+ , +—(1 + 8663 /1 1 + 8663 /1^ + ^—^ — + C'(/i^), (4.105 

80h 64 ^ 288 ^ ■i' ^ 4503 ^ v ; > v j 

where constant 63 is fixed by the condition g3{ho) = 0. Using this in p = z/h + z^gs/Rf and 
evaluating at z = 6, p = H and h = hic{5), we find 

Note that we have only kept the terms of order l/i?f in the above expansion. Now for small 
h, the integrands in /j's behave as 



f 5 + 6/^2 + 6^2) o 1 



40/i/i2v^l + /i2 + /j2 20/i2 20 

g3/^ ^ 63 1 + 140&3 (1 + 3663) (1 + 3663) 3 g.,5^ 

/i3 8O/1 64 288 V J 



~ h^Klh^ + 0{h^). (4.107) 



1 + /i2 + (^/ij 2 
/i/iVl + h? + h? 
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Using these, we can make 6 integrals in Ii and I2 finite by separating the divergences: 
h = 




^ dz / 1 1 



^4 



hic{z) ho 




(4.108) 

1 



40/iVl + + /i2 20/i2j 20 J^^ z^ \hiciz) ho 
= l's + ll (4.109) 

In (4.108), I[ and I2 are first and second integrals in Ii and similarly in (4.109), I'^ and I'^ 
represent first and second integrals in 12- Now we can further differentiate /j's with respect 
to 5 and Taylor expand the expressions to separate the divergences 



dl[ _ 1 r\ Jl + h^ + h? Kl{%^H^ + Rl) 



dV, ^ g A _ _ ^ _ (1 + I4O63) _ 1 + 3663 ^.^0. 

dS 55 V Rj J hoS^ 80i?253 64HR15 ^ "^^^ ^ ' ^^■'''> 



' - I V . . + ' ^ + , (4.112) 



dS 6^ Jh, \^40/iVl + /i2 + A2 y 20^9 

d/^ _ 3/^ / 63^/2 \ 3 3(3 + 14063)i^ 



- 20^3 ;+^+ \600i^2/ (4.113) 

Also, we can simplify the boundary terms in a similar fashion and find that 
dh _ ^hH^ \ hH'^ ^ 463^^2 ^ ^^(1 + I4O63) 



dd \ Rl ) (55 V Ri J 8O53 

+ i - (1 + 48063 + 2680061) ff \ 1 ^ ^..o. .4 i^^. 

Using (4.100) and (4. 108)- (4. 115), we find new divergence of order 1/5^ and 1/5, and the 
universal term in the entanglement entropy 

S^'^^.R^.s^ = '-^^lf^{'^ + ^{jll)^og{6). (4.116) 

Once again, we find that leading contribution in logarithmic divergence disappears, which is 
consistent with the fact that there is no contribution from the singularity with odd dimensional 
locus. In the next subsection, we will consider some more geometries with conical singularities 
to push our hypothesis. 
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4.2.2 Crease c„ x S 

In this section, we will calculate EE for conical singularities of the form c„ x S"". We will 
mainly consider following singular geometries: ci x S^, ci x S'^, ci x S"^, C2 x and C2 x 5^. 
The case with {n,m} = {1, 1} will turn out to be trivial and it will be straightforward to 
see that there is no new contribution to the log^(5 or log divergence. For {n,m} = {1,2}, 
we will find that there is a log^6 contribution in EE which had disappeared when the locus 
was taken to be flat in (4.29). For {n, m} = {1,3}, we will see that there is no log^d 
contribution from the singularity as the locus is odd dimensional. Finally for {n, m} = {2, 1} 
and {2,2}, we will find that entanglement entropy contains a logarithmic divergence. In 
{n,m} = {2, 1}, this contribution is actually coming from the smooth part of the surface as 
d = 6 and logarithmic contribution from the trace anomaly is non-zero. In {n,m} = {2,2}, 
the logarithmic contribution comes from the singularity and this is consistent with the idea 
that for even dimensional curved locus, the singularity in odd d will contribute through a 
logarithmic divergence. 

To begin with, we consider the simplest case with m = 1. In this case, the background 
geometry for the CFT is x S^. Then, the dual bulk geometry is given by 

ds^ = — {dz^ + /i(z) (di^ + dp2 ^ p2^Q2 ^ p2 sin2(^)(i02^ + f2{z)Rl df) , (4.117) 

where /i = 1 + 0{1/ Rf) and /2 = 1 + 0{1/Rf). For this bulk, we consider the singular surface 
given by G [0, il], (/) G [0, 2ir], p G [0, H] and G [0, 2tt]. For the minimal area surface which 
gives us the entanglement entropy, we assume that p = p(z, 9). Then, EE is given by 

55 I 



^4 y '^^'^^ ^4 VP +P + (4.118) 

where p = dgp and p' = dzP- Using this, we can find the equation of motion for p{z, 9) to be 

= 2zhh sm(^)p2 (p2 + p2) p" + 2z/2 sin(0)p2 (i + f^p'^) p _ 4^/1/2 s\n{9)p' pp' p' 
+2z/i/2/>2 {-2s\n{9)p + cos{9)p) p'^ + h sm{9)p^ {3zhf[ + fi (-8/2 + zf^)) p'^ 
+2zf2 {cos{9)p (p2 + p2) _ sm{9)p [2p^ + 'if?)) 

+ smi9)p^ [Azhf[ + /i (-8/2 + zf^)) {p" + (?) p' . (4.119) 

As the conical singularity has a one dimensional locus, we expect that there will be no 
logarithmic contribution from the singularity. Now we substitute p = po + pi/Rf in (4.118) 
and find the leading order correction to the entanglement entropy. However, when we use the 
equation of motion for pQ and simplify the entropy functional, the contribution depending on 
the coefficients of leading order terms of /i and /2 will vanish because /i = 1 + 0{1/Ri) and 
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f2 = 1 + 0{1/ Rf). So we find that simplified expression has following form 



cix5i 



dz dO 



Popopi 



sm{9)po / 1 d / 

-^^JPo + Pod + Po ) + I ^v^fTpfTITpp) 



+- 



1 d 



Pop'oPi 



(4.120) 



Now we can insert the ansatz pQ = z/h{6) and pi = z''^g^{6) in the functional and find that 

53 sin(6')A 



Icix5i 



^ dz /■^i-(^) dh 



z^ 



sm 



1 + /l2 + /l2 
7l3 



1 



/llc(z) 



+ 



/"''''^^^ 53sin(0) 

,n Ao /i /l2x/l + /i2 + /j2 



(4.121) 



Also from (4.119), the equations of motion for h and (73 reduces to 

= /i/i (1 + /i^) sin(e) + cos{6)hh^ + (4 + /i^) sin(e)/i2 

+ cos(6') (1 + /i^) /lA + 2 (2 + 3/1^ + /i^) sin(6l) 
= /i2 (1 + /i2) 2 sin(6')^3 + /i (1 + /i^) (2 (8 + 3/1^) sin(6')A + cos 

'2(1 + /i^) (10 + 8/1^ + /i^) sm{9) + 3 (8 + 3/1^ + /i^) sin(6')/i2 + 2 cos(6l)/i (4 + /i^) ) 93 



(4.122) 



(0)/i(l + /i2 + 3A^))53 (4.123) 



Now, we find that gs is not sourced by h and hence, a homogeneous solution gs = will be the 
exact solution for this case. This implies that excluding the first term in the above expression, 
all the other terms will vanish. As g^ is zero, there will not be any new contribution to the 
limits of the integrations either. This result is consistent with the idea that singularity will 
contribute in EE only if the locus is curved and even dimensional. 

As a next example, we consider the singular geometry ci x 52 in CFT background R^xS'^. 
For this case, the bulk metric is given by 



ds' 



L2 



{dz^ + fi{z) {dt^ + dp^ + p^de^ + p2 sm^{e)d(l)^) + f2{z)Rl dnj) , (4.124) 



where dOj = d^o + sin2(^o)'^'?o metric over unit two-sphere and 



/i 

/2 



1 + 



1 



20Rl 



+ 



WORf 
7z' 



+ ... , 



5Rl SOORf 
Now the entanglement entropy is given by 



+ ... . 



lcixS2 



Wtt'-^L^RI 

4 



dzde 



hh sin(0)pVp2 + p2 (1 + 



(4.125) 



(4.126) 
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and equation of motion for p{z, 6) becomes 

= 2z/i/2 sin(e)p2 (p2 + + 2z/2 sin(0)p2 (1 + ^ _ 4 ^ sin(e)pVp' 

+2zcos{9)f2P + p2 (1 + ^ sin(0)p( - 2 2/2(2^2 + 3p2) 

+2 p (2 z hf[ + /i(-5 /2 + z /^)) (p2 + p2)^/ _ 4 ^ 



(3 z f2f[ + 2 A (-5 /2 + z /0) 



/3 



(4.127) 



Now we can plug in the ansatz p = pQ + pi/Rf in the entropy functional and equation of 
motion. We can Tayfor expand the entropy functional in Ri and keep the leading correction. 
In this leading correction, we can use the integration by parts and the equation of motion for 
po to simplify the entropy functional to 



Sr 



cix52 



dz 



dO 



d 



sin((9)P(JPoPi 



+ 



levr^L^ 



dz Vz5iiVp2 + ^2(i + ^,2) 

We can further use po = z/h{0) and pi = z'^g^{9) and find that 

IQtt'^L^RI 



dz- 



sin(g)po {QpI + pI (6 + 5p[,=^)) 

mRlz^M+M^^W) 
s\u{e)popopi 



(4.128) 

=0-e 



6|cix52 



^1 



+ i^f + i?2 



where 



s /■''i-(^) sin(0)\/l + /i2 + /i2 

<5 ^feic(^) sin(6l) ('s + 6h^ + 6h^ 

dh — =■ 

ho 40 A /i Vl + /i^ + ^2 



dz gssm{6)h 



h=h-ic{z) 



(4.129) 

(4.130) 

(4.131) 
(4.132) 



Note that term with derivative with respect to z in (4.128) vanishes. Now using the above 
ansatz in (4.127), we can also find the equations of motion for h and 53. For these equations, 
we can make a change of variable from 6 to y = sin(0) and find 



2/(1 - y')h{l + h')h + (1 - y^)\h' + 2/(1 - y'){5 + h')h 



,2x; 2 



(4.133) 



y{l-y')h'{l + h'yg3 + h 



+(1 - 2y^)h{l + h^)h + 2/(1 + h^){5 + 2h^) 
62/ (-1 + 2/^) (3 + 4/i2 + /i4) A 



-l + 2y^)h(l + h^f 



3 (-1 + 2/2)^ /i (1 + /i2) /i2 + ( _ y (1 + /j2) (25 + 21h^ + 2/1^) (4.134) 



-3y (-1 + y2) (9 + 2h^ + /i^) + 2 (-1 + 2/^)' /i (4 + /i^) ^3)53=51, 
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where Si is the source terms and it is given by 

^ _y{l + h') (10 + 9h') + 4y (1 - y') (l + 2h^) + {l - y')' hh^ 

- 20h • ^^■^^^> 

Now to separate the logarithmic divergence, we want to find the asymptotic behavior of inte- 
grand in terms of h, where h ^ 0. So we invert equation (4.133) using h{y) = — ((l — y^) y" + 
yy''^)/y'^ and h{y) = \/l — y^/y' , where on the right hand side we have y = y{h) and 
y' = dy/dh. We can also change the independent variable in (4.70) from y to h. Apart from 
previous two relations, we also use 

Finally, we can rewrite (4.133) and (4.134) as 

h{l + h') y (1 - 2/2) y" - (5 + 7h^ + 2h^) yy'^ (4.137) 
-h{l + h^) {l-2y^)y'^-{5 + h')y{l-y')y'-h{l-y')' = 
(1 + h'')\{l- y'g; + h{l + h^) (2h (1 -y^f+ (13 + hh") y (l - y^) y' 

- (1 + h^) (5 + 2h^) yy'^^gs + (2h (4 + h^) {l - y^)^ + 3 {9 + 2h^ + h^) y {l - y^) y' (4.138) 

- (1 + h^) (25 + 21h^ + 2h^) yy'^'^gs = S2 , 

where now = dgs / dh and 

/, (1 - _ 4 (1 + 2h^) y (1 - y2) y'-{l + h^) (lO + 9^^) 

= 20h ■ (^-^^^^ 

Now we can try to solve these equations perturbatively in terms of h near the boundary, 
where h is small. As y = sin(r2) at /i = 0, we find the solution 

y = sin(O) - ^os(l^)cot(0) ^2 ^ ±_ (2csc(J^) - 7csc3(^^) + 5sin(fi)) /i^ + . . . . (4.140) 
8 512 

Using this in (4.138) and solving it perturbatively, we find that 

g^ = -^+^i±^^^^h log (h) + b3h+... , (4.141) 

y.i 20/1 ^ 96 s V ; -r d -r , V ) 

where 63 is a constant and it is fixed by the condition 33(0) = 0. Now, we evaluate the 
expression p = z/h + z^g^/Rf at z = 6, p = H and h{z) = hic- In this relation, we use 
(4.141) and invert it to find 
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Now, we can use (4.140) and (4.141) in integrands of Ii and I2 to find their behavior near 
the boundary: 



sin(e)Vl + /i2 + A2 
hh^ 



sin(O) 3cos(17)cot(J]) 3/i(13 - 19 cos(2J7)) cot2(17) csc(17) 
P~ ^ m 4096 



3sm(n) cos(17) cot(J7) /i(67- 157cos(2J7))cot2(mcsc(J7) ,^,,3^, 
H ^Ti ^ttt:^:; 1- CS(M^) 



40M3^1 + /,2 + /,2 20/i3 ' 64/1 81920 

Using these, we can make the integrations in Ii and I2 finite and write in the form 



ho 



sin(6l)\/l + /i2 + A2 sin(O) 3 cos(^^) cot(Jl) 
3/1(13 - 19cos(2Q))cot2(Q)csc(Q) 



+ 



4096 



dz f sin(r2) 

73 



I'i + I'2 

ho 



1 1 \ 3 

+ j^j + ^ cos{n) cot{n) iog(/iic//io) 



+ 



3 {hi - hl^) (13 - 19cos(217)) cot2(Jl) csc{n) 



8192 



(4.144) 



dh 



sm{9) (5 + 6/i2 + 6/i2 j 3 ^^g(j^) ^Q^(j^) 



+ 



+ 



40i?2A/i3^1 + /l2 + /j2 20/i3 64/i 

(67 - 157cos(217)) cot2((7) csc(17)/i" 
81920 

c/z /3sin(^^) / 1 1 



40 



cos(J7) cot(O) ^ , 
hi hlJ^ 64 
(^0 + ^ic) (67 - 157cos(217)) cot2(J7) csc(Q) 
163840 



(4.145) 



In (4.144) and (4.145), ![ and I2 are first and second integrals in Ii and 13 and I'^ are first 
and second integrals in I2. Now we can take derivatives of /''s with respect to (5, and then 
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Taylor expand the terms to find: 



dl[ _ 1 

~d5 ~ P 





dh 

ho 



sin(6')\/l + /i2 + /i2 sin(n) 3 cos(Q) cot(n) 3/i(13 - 19cos(2Q)) cot2(Q) csc(f^)" 



+ — — + 



/i/i3 /i3 32/i 4096 

3 (csc(J^) + 11 csc^{n) + 3 csc^{n) - 15 sin(f])) 



8192if4 

dl!^ _ H'^sm{n) csc(O) + sin(J)) log(5/i7) 

~ 26^ ^ 96^2 ^ 



-6\og{6/H) + 



d5 6 





dh 

ho 



+ 0{l/6), 



sin(g) [5 + 6/i2 + 6/i2J 3sin(^) cos(l]) cot(17) 
(67 - 157cos(21^)) cot2(J7) csc(l^)/i 



+ - 



81920 

< _ 3g^sin(^]) cos(O) cot(l]) log((5) , 

d5 4053 64 6 +^^^/'^^' 

^ H^smjn) ^ (3-cos(20))csc(17) log((^) ^ ^ _ 
(i5 20(^3 192 (5 

Now using (4.144)-(4.146) in (4.129), we find that the double log contribution in the EE 
bGConiGs 

Note that there are other contribution to the double log term but they are of higher order 
in As the only other dimensionful quantity in the problem is the IR cut-off H, these 
terms will be of the form 0{H'^ / B^). Interestingly, all such terms, which scale with H are 
contributions from the smooth part of the entangling surface. Hence (4.147) is the complete 
contribution from the singularity alone. In section 4.1.2, we saw that there was no such double 
logarithmic term when the locus of the singularity was either flat or it is odd dimensional. 
Hence, (4.147) is also consistent with the idea that similar to (2.6), generically the contribution 
from the singularity should be of the following form 

Sumv ~ / d^-^yVhlTV^] \og{5f. (4.148) 

Here a is the 2m-dimensional locus of the singularity and [TZ^] is the curvature invariants 
with m powers of the curvatures. 

Having seen the appearance of the log^J divergence for even dimensional curved locus, 
now we turn towards odd dimensional locus. So we consider EE for geometry ci x 5*3. For 
this case, the calculations proceed in a similar fashion and we find that, near the boundary 
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y = sin(^) and 53 in terms of h are given by 

y(h) = sm(n) cos(ri)cot(17)/i - — — — — 

y\ J ^ ^ 10 ^ ^ ^ ^ 6000 

(373csc(J7) - 1853csc3(17) - 889csc5(17) + 2369sin(0)) „ 

450000 ^ 

= - - "'^^-"^to"""'^'"' + + °(*'' • (^■"») 

where 63 is a constant and it is fixed by the condition that 53 has a minimum at ^ = 0. 
However, we will see that this constant is zero. If we calculate EE, we find that there is no 
log^S term but there is a log term in EE at the order next to the leading order in Ri: 

Q „4 r 6 o 

Now, we are going to argue that this logarithmic term is not coming from the singularity. We 
can see that, if we set = vr/2, the conical singularity in the entangling surface disappears. 
However, the logarithmic term in (4.150) is still non-zero. Hence, this contribution is actually 
coming from the non-singular part of the entangling surface. Further, as we are considering 
the EE in odd dimensional CFTs, there should not be any log term from the smooth part 
of the entangling surface. Hence, the contribution (4.150) should be zero and this is only if 
63 = 0. This result further assures that a singularity with odd dimensional locus does not 
contribute through a log or log^S term. 

We can further calculate EE for singularity C2 x in d = 6. In this case, the singularity 
has an odd dimensional, curved locus which does not contribute in the universal term. How- 
ever, as d is even, we find that the smooth part of the surface contribute through a log term 
and it is given by 

(4.151) 

Further, there is a new divergent term of the order 1/5 in this case. Note that this term 
matches with the logarithmic term in EE for C2 x in (4.30). Finally, we give results for 
the case where a conical singularity in odd dimensions have a curved, even dimensional locus. 
In background R^ x S'^, we can consider the singular geometry C2 x which contains a new 
divergence of order 1/6^. Further, in odd d the smooth part of the entangling geometry does 
not contribute through a logarithmic term. However, a singularity will contribute through a 
log((5) term if the locus of the singularity is even dimensional and curved. Precisely this is 
what we see in EE for C2 x 5^ and find following universal term in EE 



clog I 32 7r^I/^ 



^« „ /Co 7sin2j7 cos^O 
dh\ ^ — + 



+ 



h 60 /i4 25 /i2 
Isin^Q cos^Jl' 
180 hi ~ 25 ho 



(4.152) 
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where h = dh/dO and 

sm2(0) (7A2 + 7/i2 + 6) 

= ^ ^ . (4.153) 

60 /i4 V/i2 + /i2 + 1 

Note that similar to previous cases, h is defined such that p = z/h{6) + g^{6)/B^ and it is 
the solution of following equation of motion 

h{l + h^)h + 2 cot2(6l) /i /i^ + 2 (/i^ + 3) /i^ + 2 coi{9)hi (l + /i?) A + 3 (2 + 3/1^ + /i"^) = , 

(4.154) 

with ho = h{0). The log divergence (4.152) is a contribution from the singularity and it is 
non-zero because the locus of the singularity is even dimensional and curved. So the examples 
in this section reaffirm that an extended singularity contributes in the cut-off independent 
terms through log or log^S terms only if its locus is even dimensional and curved. 



5. Universal terms and the central charges 

In previous sections, we calculated the EE for various surfaces and found that singularity pro- 
duces new log 6 and log^d terms in the EE. As it has been seen that the regulator independent 
coefficients contain central charges when considering smooth surfaces, it is natural to ponder 
what function of central charges appears in these new universal contributions arising from the 
singularities. In the calculations above, we have been working with CFTs which are dual to 
the Einstein gravity. For these CFTs, all the central charges are equal and there is no way to 
distinguish these in the universal term of the EE. It has long been known that to construct a 
holographic model where the various central charges are distinct from one another, the grav- 
ity action must include higher curvature interactions [32]. In part, this motivated the recent 
holographic studies of Gauss-Bonnet gravity [28] — for example, see [30, 31]. Hence we apply 
this approach here as a first step towards determining the dependence of the new universal 
contributions on the central charges of the dual CFT. In particular, we will calculate EE for 
some singular geometries in Gauss-Bonnet gravities below. 

5.1 Singular embedding 

In this section, we will discuss EE for cone geometry in d = 4, 5, 6 dimensional CFTs. We 
will first calculate the EE for d = 4 systematically and discuss the results for other cases. 
For the Gauss-Bonnet gravity in d -|- 1 dimensions, the action is given by 



p I d{d-l) XL^ 
,^+^^+(d-2)(d-3)'^^ 



(5.1) 



where 

Af4 = Rabcd.R"'^'^'^ — ^RabR""^ + (5-2) 

corresponds to the Euler density on a four-dimensional manifold (and hence we must have 
d > 4 here). We have introduced L as a canonical scale in the curvature-squared interaction 
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so that strength of this term is controhed by A, a dimensionless couphng constant. We may 
write metric for the AdS vacuum solution as 



f 2 

Here the AdS curvature scale L is related to canonical scale L by the following relation 

L2 = l2//^ where = ^"^^^"^^ . (5.4) 

Now in Gauss-Bonnet gravity, the holographic EE is again determined by a minimization 
problem over surfaces m in the bulk geometry which match the entangling surface S at the 
asymptotic AdS(i+i boundary. However, the functional that is minimized is now given by 
[22, 27] 



2 



2AL _ 

1 -I TZ 

{d-2)id-3) 



(5.5) 



where TZ is the Ricci scalar for the induced metric h. 

The bulk theory can now be characterized in terms of two dimensionless couplings, the 
ratio L/ip and the new coupling constant A. As a result, there are two distinct central charges 
which characterize the boundary CFT dual to Gauss-Bonnet gravity. A convenient choice for 
these which applies for any value of d is, e.g., see [30]: 

Ct = T^^TT^ If) [1 - 2A/00] , (5.6) 
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The physical role of these charges in the boundary theory is as follows: The central charge Ct 
controls the leading singularity of the two-point function of the stress tensor and appears 
as the universal coefficient in the entanglement entropy of a sphere S"^"^ [33^ 34], xhe latter 
has also been shown to satisfy a holographic c-theorem in arbitrary dimensions. Below we 
will calculate the EE for various singular surfaces and would like to see if the universal terms 
have some simple dependence on these central charges. If we consider the case of d = 4, the 
above expressions simplify to 

c = Ct = 7r2 — (I-2A/00) and a = = vr^ — (1 - 6A/00) • (5.8) 

Here, c and a are the standard central charges that appear in the trace anomaly or entangle- 
ment entropy, e.g., eq. (1.5), of the four-dimensional boundary CFT. 

Now for calculation of the holographic EE, the cone geometry is defined by p G [0,-ff], 
Q G [0,r2] and <^ G [0,27r]. We consider the induced coordinates to be {p,6,(p) and radial 
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coordinate z = z{p,9). For this case, the induced metric becomes 



h 



(1 + (dpz)') 
^dpzdgz 




^dpzdgz 

{p' + (dez) 









(5.9) 



For this metric, the expression for Ricci scalar TZ in (5.5) contains the terms like dpZ and d^z. 
However, it is straightforward to see that the equation of motion is still second order. This 
is because Gauss-Bonnet term is topological in nature. Further, we impose the UV cutoff at 
z = 5 and define e{p) such that at 9 = 0, — e, z{p, — e) = 5. As the background geometry 
has scaling symmetry and apart from p, there are no dimensionful quantities in the problem, 
the solution for z should be of the following form 



ph{9) 



(5.10) 



Here h{9) is a function such that h{Q) = and /i(0) = 0. Also, the maximum value of h{9) is 
/i(0) = /iQ. By plugging this ansatz in equation of motion for z{p, 9), which we get by applying 
the variational principle on entropy functional (5.5) with d = 4, the equation of motion for 
h{9) turns out to be 



= /i(l + /i2) (^(l + /i2)sin(^)(l + 4A/oo)+6Acos(0)/i/oo/i + sin(e)(l-2A/oo)/i^j/i 
+ cos{9)h (1 - 2A/oo) /i^ - (3 + h"^) sm{9) (-1 + 2A/oo) 

+2 cos{9)h (1 + + A (1 - 2/i2) /^) - 3 (1 + /i^) sm{9) (-2 - /i^ + 2A (l + h^) foo) h? 
+ cos{9)h (1 + h^) 2 (1 + 4A/oo) A + (1 + 2 s\u{9) (3 + 2/i2 (1 + A/oo)) . (5.11) 

Further, we can simplify the entropy functional using this equation of motion and find that 



ICl 



H 



dp 



5/ho P Jho 



5/p 



dh 



sm{9)Ci 

Co 



(5.12) 



where we have changed the integration from 9 to h and 

= (l + h^ + ( sin(^) (l + h^ + h^) - 2A/oo ( {h^ cos{9) cot{9) + 4 sm{9)) 



+2h^ cos{9)h + (2 + h^) 2 sin(6') j j + AX^ {h^ cos{9) cot{9) + 3 sin(6')) 
-2h (3 - h^) cos{9)h-^ + /i^ (2 + + 2cos(20)) csc{9)h^ + 2h^ (l + h^) cos{9)h 
(1 + /i2) sm{9) 



(5.13) 



k^hy 1 + /i2 + /i2 sin((9) 1 + /i^ + + 2A/oo 3h cos{9)h + sin(6') (2 + 2h 



Now we want to make h integrand in (5.12) finite. For that, we define y = sm{9) and 
find y and h in terms of h near the asymptotic boundary. For that, we convert (5.11) into 
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the equation of motion for y with independent variable h. Solving this equation of motion 
perturbatively, we find that near the boundary 

y{h) = s'm{Q) — -h^ cos{n) cot(O) + — /i^ log(/i)(3 - cos(21^)) cot^(17) csc(n) + co/i^ + .(5..],4) 
4 64 

where we have used the condition y(0) = sin(O) and cq is a constant which is fixed by the 
condition that y'(/io) = 0. We can now assume y = y{9) and h = h{9) in (5.14) and then 
invert it to find 
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Now using (5.14) and (5.15), we can find that near the boundary 

sin(0)A _ sin(^7) (1 - 6A/00) ^ cos{^) cot(17) (1 - 2\f^) ^ ^^^^ 
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which can be used to write the EE as 
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In the last term, the h integration is finite in the limit 5 — )• and hence the leading order 
divergence is logarithmic. Hence, we write 
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Now we can compare the coefficient of the log^6 divergence with the central charge (5.6) and 
find that 
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where the second equation follow from (5.8). So we find that for EE of the cone ci in d = 4, 
the new universal term arising from the singularity is directly proportional to the central 
charge Ct and has a relatively simple dependence on the opening angle il. We will return to 
discussing this result in more detail in section 6. 

We can also calculate the EE for five-dimensional cone in Gauss-Bonnet gravity. For this 
case, the action for bulk geometry is (5.1) with d = 5 and canonical scale L is related to L by 
relation (5.4). Now the complete expression of entanglement entropy is given by (B.8) and 
the universal term is 

,,iog I 8 TT^Z^ / (2cos^(17) (2 - 7A/oo) hi -3sm\n) (1 - AXf^)) f^o r\. 

s,\ = -^ + i '^ZI j 

where C3 and £4^ are given by (B.9). In this case, we can further compare the universal term 
with the central charges (5.6) and (5.7). However, we observe from expression of £3 and £4 
in (B.9) that there are terms of order ©(A^) in the universal term. This implies that the 
expression (5.20) is not a simple function of the central charges Ct and 05 and in particular, 
it is not a linear function. 

Further, we calculate the EE for cone in d = 6 and find that the regulator independent 
term is log^d and it is given by 

.log^i _ 12 tt^LS cos(O) cot(^7) f93 - 190A/oo 3 - 2A/oo \ . 



6 lc3 £5 \ 8192 8192 V ; / & V / A ; 



We have given the complete expression of EE in the appendix B in eqn. (B.IO). In this 
expression, we would like to compare the $7 independent coefficients with the central charges 
of the GET. So there are two terms, (93 — 190A/oo) and (3 — 2A/oo), which we can express in 
terms of Ct and Og. Using equations (5.6) and (5.7), we find that the universal terms in EE 
for C3 can be written as 



C,log2 



3 cos(r2) cot(r2) 
^=3 ~ 1024 



90 Ct + 3 - 6 Ct - 3 og cos(20) log^{6/H) . (5.22) 



Hence we see here that in d = 6, the new universal term for C3 depends linearly on both of 
the central charges, Ct and Og, but it still has a relatively simple dependence on the opening 
angle Q. 

5.2 Singularity with a curved locus 

In this section, we will repeat the calculation of EE for geometry k x S'^ in d = 5 GET, which 
is dual to the Gauss-Bonnet gravity. For this case, the bulk action and EE are given by (5.1) 
and (5.5). We further consider the metric ansatz 

= _ (^^2 ^ ^^^^^ (^^^2 ^ ^^2 ^ ^2 gij^2(^^^ ^ Rlf2{z)dnl) , (5.23) 
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where dQ,2 is line element over the two-sphere. Further, similar to (5.4), = L'^ / foo where 
/oo = (1 ~ \/l — 4A)/2A. Here /i and /2 are functions of z and one can use the Fefferman- 
Grahm expansion to find their values near the boundary. 

Now, to calculate EE, we first use Fefferman-Grahm expansion and find that 

= 1 + ^+ ^^(51-58A/^) 
' ^ 12/^2 ^ 1728i?4(i - 2A/oo) ' 
f(\ 1 z^(15-2A/oo) 

= '-4i?f- 576i?t(l-2A/..)+---- ^'-''^ 

Once again, we choose the parametrization p = p{z, 9) and similar to (4.35), we can find the 
induced metric for the entropy functional. Here, we restrain to give the complete expression 
for the entropy functional but the important point to note is that now there are terms with 
higher derivatives, like p", p and p' in the entropy functional. However, in spite of this, the 
equation of motion for p[z, 6) is still second order and that is because of the topological nature 
of the Gauss-Bonnet terms. Now we can simplify the entropy function by inserting the ansatz 
P = Po + Pi/R\ and then series expending it for large Ri. We find that up to the leading 
order, the EE is given by 

•S'slfcxsa = Ji / / (^^i^o{Po) + -j^y^N{po)+Cf{po,Ci,C2)+p'iCi{po) 

+fhC2ipo) + Pi^siPo) + PiA(po) + /5iA(/5o) + Pi^eiPo)) I , (5.25) 



where Cq{pq) is the term which comes from the limit i?i — )• oo in the original lagrangian. The 
term Cm{pq) is a new term which doesn't appear for A = case. This term comes from the 
contribution of sphere in the Ricci scalar in (5.5). The term Cf{pQ, ci, C2) is the term which 
is independent of pi and comes from the leading order corrections in /i = 1 + Ciz'^/Rf and 
/2 = 1 + C2z'^/Rf, where ci = 1/12 and C2 = —1/4. Further, terms with Ci{pq)-jCq{pq) are 
independent of ci and C2 and are linear in pi, as it is written. Now, we can write 

p[Ci = dziCipi - C[pi) + C'Ipi , 
p[C2 = dz{C2Pi) - de{C'2Pi) + Ci pi , 
P1C3 = de{C3pi) - deiCspi) + Cspi , 
piCi = dz{Cipi) - C'^pi , 

P1C5 = d0{C^pi) - C'^pi , (5.26) 
where prime and upper dot denote the partial derivative with respect to z and 9. Using these 
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in (5.25), we write 

P J Zfn "/ \ 1 

+Pi{C'l + £2' + £3 - ^4 - A) + dz{£ipi - C'lPi + C2P1 + A/Oi))^ 

+ dz (^-C2Pi + C-ipi - C3P1 + C5pi^^_^^^ , (5.27) 

where the coefficient of the pi vanishes from the equation of motion of po and we have 
performed the integration over 6 in terms in the last hne. Note that we can not integrate 
over z in any term as e = e{z) and both 9 and z integrations don't commute. Now we further 
insert the ansatz po = z/h{6) and pi = z^g^iO) in the above entropy functional and in the 
equations of motion for po and pi to find the equations of motion for h and g^,. Once again, 
we restrain to give the complete expression of the equations of motion as they are not very 
illuminating. Similar to (4.45), at this stage we find that the terms with the partial derivative 
with respect to z vanish. Finally, the entropy functional reduces to 

„, leiU'flf /' dh ( Co(li} ^ I CN{h) + C,(h,c,,C2)\ 

s^Ls^ = X [— + M ; j 

(5.28) 

where we have defined Cq = Co{h)/z'^, Cjy = Cj\[{h)/z, Cf = Cf{h,ci,C2)/ z and expressions 
for £'s are given by (B.13) in appendix B. Now we can solve the equations of motion for h 
and gs near the asymptotic boundary to find 

. _ _ai _ 2ai _ (1 - 6/ooA) _ /^^(I-IO/^A) 

+ 2ai (1 - 2/ooA) 2ai(l-2/ooA) ' ^^"^^^ 

^ b3 9 -26/00A + 79263 (1-2/00 A) 4/i(3-ll/ooA + 54fe3(l-2/ooA)) 

504/1 (1-2/00 A) 567(l-2/ooA) ' 

where oi and 63 are constants which are fixed by ensuring that both h and 53 have extrema 
at 6 = 0. Note that here ai is related to a quantity which is conserved along the 9 translation 
similar to (4.46). Using these solutions near the boundary, we first find the value of hic{z) at 
the UV cut-off z = 5. To do that, we use the above solutions in the ansatz p = z/h + z^g-^/R^ 
and inverting the relations iteratively, we find 



hiciS) 



1 + hE\ s + (9-26/ooA + 792fc3(l-2/ooA)),5 
H Rf J 504ifi?f (1 - 2/ooA) 

4(3-ll/ooA + 5463(l-2/ooA))^5 

567H3RI (1 - 2/ooA) 



+ . . . . (5.30) 
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Note that above relation reduces to (4.49) for A = 0. Now we use (5.29) to study the behavior 
of integrands in (5.28) near the asymptotic boundary, where we have /i — >• : 

^0 1 - 4/ooA ^ ^^^4) ^ 



h h? 
CN + Cf 15-68/ooA 



+ 0{h^). (5.31) 



h 72/i2 

Using these, similar to (4.52) and (4.54), we can make the integrands in (B.13) finite. We 
can break the terms in following components 

[' dz ['"'('> (it, 1 - if^\ \ 



and then take a derivative with respect to 5 and find following series expansion: 
dh _ bsH^\ H{l-Af^X) 1-4/ooA 



H{1- 4/ooA) (9 - 26/ooA + 79263 (1 - 2/ooA)) 



dh _ I f\,,( CN + Cf 15-68/ooA 



504 (1-2/00 A) RfS^ 

^ - _ /^i _ g (15-68/ooA) 15 - 68/ooA ^ 

dS V Rl J 7252 + j2ho5 ^ 

d CB\e=n-e ^ A 3b-iH^\bsH^l-4f^X) 



ds J, z V ^? y 

/ _ 463^\ iJ (1 - 4/ooA) (9 - 26/ooA + 792^3 (1 - 2/ooA)) ^ 
R\ ) 504(l-2/ooA)52 

Note that we have done the same with the boundary term in (5.28) too. Using these relations, 
we can read off the logarithmic term in the EE, which is given by 

log, 16^2^4 { f\^( CN + Cf 15-68/ooA ^ , 15-68/ooA V 



(5.37) 
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where and Cf are given by (B.13). Now we can try to compare this cut-off independent 
term with the centrai charges (5.6) and (5.7). However, by looking at the expression of Cn 
and Cf, we find that there are terms of order O(A^). Hence the universal term is not a simple 
(and particularly linear) function of the central charges. 

Now, it would have been interesting to further investigate the universal term for the 
geometry ci x S'^. In this case, the singularity has a even dimensional, curved locus and we 
will get a log^(5 divergence. So it would have been very interesting to see if this contribution 
from the locus is a specific central charge. However, this calculate is tedious and demands 
more patience and ingenuous simplifications. 

6. Discussion 

In this paper, we have used holography to study entanglement entropy for various singular 
surfaces in higher dimensions and our results are summarized in table 1. In particular, in 
section 3.1, we considered cones in various dimensions. This analysis suggests that for CFT's 
in an odd number of spacetime dimensions, an additional universal contribution appears: 

= Qdin) \og{6/L) , (6.1) 

where is the opening angle of the cone, as shown in figure 1. In particular, for C2 = x S"^ 
in d = 5, eq. (3.25) provides an expression for the coefficient q^{^) in the holographic CFT 
dual to Einstein gravity. We showed that this coefficient satisfies the required properties: (a) 
QdiS^ = 7r/2) = since the entangling surface is actually a fiat plane for this angle and (b) 
grf(r2) = qd(7r — Q) since the entanglement entropy of the CFT ground state is identical for 
the density matrix describing the degrees of freedom inside or outside of the cone. Further 
we found that qd{^ — )■ 0) oc 1/17 as in eq. (3.27). Again our expectation is that this behaviour 
extends generally to cones for odd dimensional theories. That is, the entanglement entropy 
acquires a universal contribution of the form given in eq. (6.1) for a cone Cd_3 = x S"^~^ 
in any odd d. 

Here we might recall that in three dimensions, strong subadditivity was used to derive 
further constraints on the form of (73 (17) [18, 19]. In particular, one shows that this function 
must satisfy q-i{^) > and q'^{^) < 0. Of course, our holographic result (3.9) satisfies both 
of these inequalities. It is noteworthy then that in five dimensions, our holographic result 
satisfies 95 (17) < and q'^{^) > 0. It would be interesting to understand if these inequalities 
are again general properties originating from strong subadditivity. 

Further we add that in eq. (3.25), the coefficient q^ is proportional to L'^/£p which can be 
identified with a central charge in the dual boundary CFT. However, with Einstein gravity 
in the bulk, all of the central charges in the five-dimensional boundary theory are identical 
and so in section 5.1, we extended the calculation to Gauss-Bonnet gravity in the bulk, 
which allows us to distinguish at least two such central charges, as described there. However, 
given the result for ^5(17) implicit in eq. (5.20) and (B.9), it appears that this coefficient is a 
complicated nonlinear function of both central charges. It may in fact be that this coefficient 
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is not determined by the central charges alone - this would then be similar to the results 
found for Renyi entropies in [30]. 

The case of even dimensions was particularly interesting. In section 3.1 with d = 4 and 
6, we found that a cone yields a universal contribution of the form: 



where il. is again the opening angle, as shown in figure 1. Again we believe this is a generic 
result for even dimensional theories. In our holographic examples, the coefficient functions 
have a relatively simple form, as shown in eqs. (3.23) and (3.28). Hence for these examples, 
it is straightforward to verify that qd{^) satisfies the same properties as described for qd{^) 
in the first paragraph above. From these results, we also see that qd{^ — )• 0) oc 1/^2, which 
parallels the behaviour seen for qs and q^. Further, in section 5.1, using Gauss-Bonnet gravity, 
we were able to verify that this coefficient is proportional to the central charge (5.6) which 
controls the leading singularity of the two-point function of the stress tensor for d = 4. In 
particular, we found the simple result (5.19) which yields 



For d = 6, we also found a simple expression (5.22) for the new universal term for C3. However, 
in this case, % depends linearly on both of the central charges, Ct and Og, and it is still a 
relatively simple function of the opening angle fi. Given the simplicity of these results, 
particularly in eq. (6.3), one might conjecture that the same form may arise for CFT's in 
general beyond our holographic calculations. 

We note here that with even d, there are also contributions proportional to a single power 
of log(5/L, however, these terms are no longer universal. Rather the corresponding coefficient 
will vary if the details of the cut-off (or the macroscopic scale L) are changed because of the 
presence of the log^ 6/L term - see further discussion below. 

Recall that the trace anomaly in an even dimensional CFT gives rise a universal con- 
tribution in the entanglement entropy with a smooth entangling surface [5, 16, 35, 33]. In 
particular, if we consider a four-dimensional CFT in a flat background, this contribution takes 
the form given in eq. (1.5). Given this explicit expression, it is interesting to compare this 
contribution for a conical entangling surface ci in d = 4 to the universal term (5.19) found in 
our holographic calculation. 

Now in d = 6, we do not know the precise expression of the contribution of the trace 
anomaly to EE, i.e., d = 6 generalization of (1.5). Hence, now we can not distinguish the con- 
tribution from the trace anomaly with the contribution from the singularity in (5.22). Hence, 
we can not say what part of the universal term comes from the singularity and particularly, 
we can not confirm if singularity contributes only through the central charge Ct- 

Of course, eq. (1.5) is only expected to apply for a smooth entangling surface and so 
cannot be applied directly to the cone ci. Hence our approach is to construct a 'regulated 
cone' ci by cutting off the cone at some p^i,^ and replacing the tip with a spherical cap of 
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radius r = p^^^ tanil, as shown in figure 4. Hence ci provides a smooth model of the desired 
conical entangling surface to which we can apply eq. (1.5) and we can consider the limit 
Pmiii — )• to recover the result for the singular surface ci. Hence working first with finite Pmin, 
it is relatively straightforward to show that eq. (1.5) yields^ 

•Suniv = cosfi iog(j/L)log(p„,i„/L) 
2 sm \ I 

+ 2a (1 -sinf]) log((5/L). (6.4) 

Unfortunately the first term will diverge if we take the limit p^in — )• 0. However, if the 
underlying CFT has been regulated with the cut-off 5, it should not be able to resolve any 
features in geometry at shorter distances. Hence it is natural to consider the limit Pmin — ^ ^ 
which yields 

5.niv = \og\5/L) + 2a (1 - sinJ7) \og{5/L) . (6.5) 

2 sm S I 

Hence with this construction, the universal contribution (1.5) contains terms proportional to 
both log^{5/L) and log{6/L). The most surprising aspect of this result is that the coefficient 
of the leading term is almost identical to the holographic result in eq. (5.19). However, 
there is a mismatch by a factor of two.^ Our interpretation of this result is that with the 
conical entangling surface ci, part of the log^(5/L) divergence should be associated with 
correlations across the smooth part of the entangling surface away from the singularity. The 
above calculation then suggests a pile up of short distance correlations in the vicinity of 
the tip. However, the full accumulation of correlations near the tip of the cone does not 
have (precisely) the form expected from the 'smooth' expression (1.5). Hence a part of this 
universal contribution should be thought of as intrinsic to the singularity at the tip of the 
cone itself. 

As already noted, if we examine our holographic entanglement entropy for ci in eq. (3.23), 
we also find a term proportional to log{6/L). The angular dependence of the corresponding 
coefficient is different from that for the log^((^/L) term and so one might envision that certain 
universal aspects can still be extracted from this coefficient. For example, in eq. (6.5), the 
second central charge a appears in the \og{6/L) term and so one might hope to extract 
this charge by studying the entanglement entropy for cones with various opening angles 0. 
However, we wish to emphasize that this log{6/L) contribution is simply not universal in 
the presence of the log^(5/L) term. For example, let us return to our regulated model of 
cone above where we argued it was natural to take the limit p^i^ — )• 5 for the radius of the 
spherical cap. The latter was motivated by observing that 5 is a short distance cut-off in the 
underlying CFT and so the latter can not resolve any geometric features involving shorter 
distance scales. However, let us note that the radius of the spherical cap was r = p^^^ tan $7 
and so even if we set p^in = 6, for small Q, the cap is effectively much smaller than the short 
distance cut-off. Hence one might instead choose 

/ S for a>w/i. 



**On the conical portion of ci, 7?. = wliile the only nontrivial component of the extrinsic curvature is 
K^^ = P sinfi cosfi. On the spherical cap, TZ — while the combination of extrinsic curvatures in eq. (1.5) 
vanishes. 

^This observation was also made in [23]. We also note that the same result appears in an alternate 
calculation of the holographic entanglement entropy Resented in appendix A. There the mismatch can be 
regarded as an anomaly arising from a singular conformal transformation. 



While such a choice leaves the log^(5/L) term unchanged, the coefficient in the log{S/L) 
contribution acquires a complicated new angular dependence. We present this discussion 
here simply to illustrate that all of the details (including the angular dependence) of the 
coefficient of the log{6/L) can be expected to depend on the precise choice of the regulator 
in the calculation of the entanglement entropy. While eq. (6.6) gives an illustrative example, 
we might observe that the same regulator in the holographic calculations is far more subtle. 
The analog of eq. (6.6) would be p^in = 5/ho in section 3.1. 



Figure 4: (Colour Online) The illustration shows log^^ in even dimensions. These results indi- 
rcgulatcd cone ci, where the tip is replaced by a cate that there is a rich variety of new geo- 



gling surfaces. However, our analysis only considered simple families of singularities and does 
not suffice to reveal the full geometric structure of these universal terms. It would, of course, 
be interesting to consider more general singularities, e.g., a crease of the form k x but 
where the opening angle varies along or where R^ was not entirely straight. Another step 
towards a clearer picture of these geometric coefficients would be to carry out the holographic 
calculations using the Fefferman-Graham expansion [36] to compute entanglement entropy 
along the lines discussed in [15]. 

In part, our motivation for these studies was the possibility that these new universal 
contributions may be used to identify the central charges of the underlying CFT. Recently 
there has been a great deal of interest in using entanglement to identify central charges 
that obey a c-theorem, in particular for odd spacetime dimensions, e.g., see [11, 23, 33, 34, 
37, 38]. In this regard, our results only point towards a clear result for even dimensions, 
namely, that the coefficient of the log^(5 contribution is proportional to a particular central 

^''Let us re-iterate that we are ignoring certain cases here, e.g., c\ x E? , where a log 5 term appears but it 
can be attributed to the smooth part of tlie entanghng surface. We also do not consider the possibility that 
the finite contributions may exhibit some new universal behaviour. 




M = 0/ 



We also extended our holographic analy- 
sis to consider creases or extended singular- 
ities in section 4. Examining the examples 
summarized in table 1, we find that the crease 
of the form k x or x i?*" creates no 
new universal contributions.^*^ However, in 
general, we found that creases can contribute 
additional universal terms, but singular lo- 
cus must have an even dimension and must 
be curved. These results suggest that these 
new universal contributions to the entangle- 
ment entropy take the form given in eq. (2.6) 
with a log 5 divergence in odd dimensions and 



spherical cap. 



metric contributions to entanglement entropy 
that can be associated with singular entan- 
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charge. Unfortunately, this central charge Ct is not the one expected to satisfy a c-theorem. 
However, our holographic result seems particularly simple and so it may be that there is a 
general derivation for any CFT, perhaps connected to the trace anomaly as for the universal 
terms identified for smooth surfaces [5, 16, 33, 35]. It is also worthwhile to investigate if 
other central charges appear if we consider, e.g., cones with a more general cross-section than 
simply S"^~^. 

In odd dimensions, our results in section 5.1 indicate that the universal terms associated 
with singularities in the entangling surface will be complicated nonlinear functions of many 
parameters in the underlying CFT. Here it must be said that we refer to these contributions 
as universal since they will be independent of the details of the UV regulator and so should 
characterize properties of the underlying CFT (or QFT more generally). However, the precise 
nature of the information contained in these terms remains to be understood. A similar result 
was found for Renyi entropies of spherical entangling surfaces in [30] . 

It would also be interesting to study these universal contributions to entanglement entropy 
in non-holographic theories. In particular, one might consider heat kernel methods for free 
field theories for simple surfaces, e.g., along the lines of [14, 40]. We have been informed 
by Brian Swingle that he found similar log^S terms in the entanglement entropy of surfaces 
with a conical singularity using a twist field calculation in the Gaussian approximation, as 
discussed in [41]. It is also of interest to extend our holographic calculations to consider 
entangling surfaces with flat faces and edges, e.g., a cube. In particular, one would like to 
confirm that the new log^J contributions persist for the 'corners' in such cases. This would 
have practical implications for lattice calculations in higher dimensions, where it is difficult 
to avoid such corners. It would also be interesting to further investigate the Renyi entropy 
for singular entangling surfaces [17]. 
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^^Setting holography aside, it was shown that for four-dimensional QFT's, there is no possible (linear) 
combination of the two central charges, c = Ct and a — a^, that can satisfy a c-theorem apart from a alone 
[39]. 

^^For example, in four dimensions, we can regard such a corner on the surface of a cube as a 'cone' with a 
triangular cross-section. 
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A. Conformal transformations and EE for q_3 

In this section, we consider an alternate approach to calculating the entanglement entropy 
associated with a conical singularity. In particular, we begin by performing a conformal 
transformation which takes R'^ to R x 5"^^^. Under such a conformal transformation, the 
conical entangling surface Cd-3, considered in sections 2 and 3, becomes simply a uniform 
cylinder RxS'^^^ in the latter background. The analogous mapping was applied to calculating 
the cusp anomaly for a Wilson with a sharp corner in N=4 SYM, e.g., see [42]. In this case, 
calculating the cusp anomaly becomes problem of determining the quark-antiquark potential 
on a three-sphere. In fact, the holographic calculation for the leading result corresponds 
to determining an extremal surface in the bulk with a 'kink' boundary condition and so it 
precisely matches the calculation of the holographic EE for a kink in d = 3 [42, 43]. 

We would like to calculate EE using the cylindrical background geometry R x S"^~^ 
and in particular, to see if we get the same log 6 and log^ 6 terms that appeared in our 
previous calculations for the cones Cd-3. The key feature in these calculations is that the 
conformal mapping takes the conical entangling surface in R'^ to a new surface R x S"^~^ in 
the cylindrical background. Now this new surface has an infinite length along the R direction 
and so this length must be regulated to properly account for the logarithmic divergences, as 
we will discuss below in section A.l when we calculate the holographic EE. We first describe 
desired conformal transformation in the CFT and then we consider how this transformation 
is implemented with a coordinate transformation in the dual bulk spacetime. 

In the flat background (2.2) with {n,m} = {d — 3,0}, we can make the coordinate 
transformations Ie = r cos ^ and p = r sin ^ and find following metric 

ds"^ = dr^ + r'^{d(,^ + sm^^{de'^ + sin^ednl_^)). (A.l) 

In these coordinates, the cone geometry discussed in section 2 translates to the surface: 
Cci-3 = {r = [0,00),,^ = 11/2,9 = $7}. We can further perform the coordinate transformation 
r = Le^l^, and make a Weyl transformation to remove the overall factor e^^l^ from the 
resulting metric. After these transformations, we find background geometry becomes R x S'^~^ 
with the metric 

ds^ = dY^ + L^{df + sin^ e {de"^ + sin^ dnl^^)) , (A.2) 

where Y £ [—00,00]. Further, the conical entangling surface above is now mapped to a 
cylinder of infinite length, i.e., 

= {Y = {-00,00),^ = 7T/2,9 = n}. (A.3) 

Next we discuss the coordinate transformation in the bulk geometry which implements 
the conformal transformation between the two boundary metrics in eqs. (A.l) and (A.2). We 
begin with the standard description of (Euclidean) AdS^+i as a hyperbola embedded in the 
following (d + 2)-dimensional Minkowski space 

ds^ = -dU^ + dV^ + {dX^f + ■■■ + (dX'^f = -dU^ + dV^ + dR^ + R^dflj^^ , (A.4) 
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where in the second expression, we introduced polar coordinates on the space spanned by the 
X*. Now the AdS geometry is defined by the hyperbola 

We can solve this constraint by writing U = \/ B? + L'^ cosh(y/L) and V = ^ E? + LP' sinh(y/L), 
in which case the induced metric on this surface becomes 

ds^ = - — ^—^dR^ + fl + ^) + W + sin^ ^ (^^^ + sin^ ^ ^^Ls)] • (A.6) 
(l + f) V LJ 

Here we have written dVl'^^_^ = d^^ + siv? ^ {dO"^ + s\v? 9 dQp^_^) . Of course, we recognize 
this metric as (Euclidean) AdS^+i in global coordinates and this geometry is dual to the 
boundary CFT with the background metric (A. 2). Alternatively, to get the bulk metric dual 
to eq. (A.l), we write R = rL/ z and U + V = l? j z and eq. (A. 5) yields U — V = z ^ j z. 
With these coordinates, the induced metric on the hyperbola becomes 

ds2 = ^{dz^ + dr^ + r^dVL\_^) , (A.7) 

which matches the 'Poincare patch' metric (2.3) with = + p^. Clearly, the boundary 
metric matches the desired form given in eq. (A.l). This bulk metric (A.7) was used in our 
calculation of the holographic EE for the cone in section 3. So now we use the metric (A. 6) 
in calculating the EE for the cylindrical surface produced by the conformal mapping. 

A.l EE for cylinders 

In this section, we will calculate the holographic EE for cylindrical entangling surface RxS'^~^ 
on the background R x S"^~^ for d = 3, 4 and 5. As we will see the calculations here are closely 
related to those already presented in section 3. We find that the coefficient of the universal 
log 6 term matches our previous results for d = 3 and 5. In d = 4, the universal term is a log^ 6 
contribution but we find the coefficient in the following does not give a precise match with 
that in eq. (5.19) from our previous calculation. In fact, the log^ S term here shows the same 
mismatch by a factor of two that we found in eq. (6.5) from considering the contribution 
of the trace anomaly on a regulated conical surface. We will show that this difference in 
the holographic results comes from choosing different UV cut-offs in the different coordinate 
systems. Although these cut-offs are natural in the particular coordinate system in which 
they are chosen, they produce different coefficients of the log^ 6 divergence. 

Now we begin with a general discussion for a d-dimensional CFT on the background 
R X S"^~^ and later focus on the specific cases with d = 3, 4 and 5. The cylindrical entangling 
surface on the boundary is given in eq. (A. 3). However, to produce a finite result, the length 
the cylinder must be regulated which we do by restricting Y G [y_ , Y^] , where y_ and y+ are 
cut-offs to be fixed below. Now the holographic EE is calculated by first finding the minimal 
area surface hanging into the bulk geometry with metric (A. 6) and if we define this surface 
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by coordinates (¥,9,0,^-3)1 we will have a radial profile R{9) = Lg{6). In particular, note 
that R{p) is independent of Y because we have translational symmetry along this direction. 
The only scale in the geometry is L, which coincides with the radius of the sphere S'^~^ on 
the boundary. With this ansatz, the holographic EE is given by 

S2d_3 / + / .Qd-3(■n^d-3. 



^d\ cylinder 



[ dY [ ' de g'^-^{sm Of-^^/g'^ + g"^ + 9^ , (A.^ 



where g = dg/dO and ^d~3 represents the area of the unit (d — 3)-sphere. Further, e is related 
to the UV cut-off in the theory. Note that the natural way to choose this UV cut-off for 
cylinder should keep e independent of Y , the coordinate along the length of the cylinder. So 
a natural UV cut-off for the minimal area surface can be R = 1? jb^ which will fix e such that 
giSl - e) = Lib. 

To compare the above expression (A.8) with our previous results, we need to relate the 
quantities for the cylinder with the quantities for the cone discussed in section 3. First we 
note that from the coordinate transformation taking us from eq. (A. 6) to eq. (A. 7), we have 
p = r and R = r L/ z = pLjz. In our previous calculation of the holographic EE in eq. (3.3) 
or (3.14), we used the ansatz z = ph{6). Hence, above relation implies that R = L/h{9). 
Comparing with our ansatz above, we see that g{9) = l/h{9) and h{^} — e) = 5/L. Now we 
turn to consider the cut-offs in the coordinate Y. First note that these are independent of 
our UV cut-off on the radial coordinate, i.e., R = L'^/6. As the conformal transformation has 
mapped the cone in R'^ to the present cylinder, y_ and 1+ should be related to the cut-offs in p 
appearing in eq. (3.14). We use the coordinate transformation Y = L\og{p/L), as given above 
eq. (A. 2), and from the range p G [S/ho, H] in eq. (3.14), we find that y_ = Llog{6/hoL) and 
y+ = L\og{H/L). Using all these results, we can now rewrite eq. (A.8) as 

2^L'^-iJld_3 dp /■'^/■^ sin'^-3(^) 



At first sight, this expression is identical to eq. (3.14) found in the previous calculation. How- 
ever, we should notice that there is crucial difference. Namely, the upper limit of integration 
over h here is a fixed constant while in eq. (3.14), it is a function of p. 

Now we consider cylinder in different dimensions one by one. In eq. (A.8), we first take 
d = 3 to consider a cylinder dual to a kink in R^ . Note that in this case 9 G [— ^2, fi] and we 
will also have an extra factor of two coming from the change in the limits of integration. We 
can further use the variable y = y^l/Zi^ — 1//iq as around eq. (3.6) and follow the steps in 
section 3. Finally, we find that for d = 3, eq. (A.9) will produce 

^""^^ \og{H/5) + ^ log(/io) - q3m log{H/6) + ... 



S: 



3\cylinder 



(A.IO) 



where ^3 is given by (3.9). We note that the universal term in EE for the cylinder here precisely 
matches with that for the kink. However, now the leading order divergence is different as it 



^^The corresponding symmetry for the cone in R'' is invariance under dilatations. 
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contains an extra factor of log{hQH/5). Of course, these divergent terms are not expected to 
be universal. 

In part, the aim of the present calculations is to show that our 'universal' terms are in- 
dependent of the details of the choice of regulator. Here we chose the natural cut-off adapted 
to the new coordinates after the conformal transformation in the boundary geometry (or a 
coordinate transformation in the bulk geometry). While the entire structure of the diver- 
gent contributions was not unchanged, our universal logarithmic term matched the previous 
calculations. Note that in eq. (A. 9), if we would have used z = 5 as the UV cut-off, using 
R = pL/z, we would have had g{Q — e) = p/6 = l/h{0, — e). In this case, eq. (A. 9) (with an 
extra factor of two) would have become 

Js/ho P J ho h?h 

precisely matching with eq. (3.3). However, choosing the new radial coordinate i? to be a 
function of p (or more clearly Y) would not have been a natural UV cut-off in this case. 
Although we have found that for d = 3, the universal term in eqs. (3.3) and (A. 10) match, 
this will not be the case in d = 4 which we discuss next. 

For a cylinder in a four-dimensional CFT on Rx S^, the EE is given by eq. (A. 9) with 
d = 4. Now following the calculations and steps in section 3.1, we find that 



S. 



4 1 cylinder 



47r^ 
^3 



L^sin^J]^ ^^^^^ L^in^n^ 

log{H/5) + log(/i, 



0) 



^ (A.12) 



cos $7 cot $7, 2/r/rTN ^^z-, /cw 

log\5/H) + 0{log{S)) 



In comparing eq. (A.12) with eq. (3.24), we see that the universal term is off by a factor 
of two. However, this log^ S contribution precisely matches with the contribution from the 
trace anomaly as given in eq. (6.5). Here we can see that the new universal contribution from 
the singularity discussed in eq. (5.19) is not invariant under the conformal transformation 
considered here. In eq. (A. 9), if we had chosen the UV cut-off on i? = pL/z = Lg{6) to be 
z = 6, the upper limit of integration would have been 5/p, instead of 5/L. In that case, the 
results for the holographic EE here would have precisely matched our previous results for the 
cone ci. 

Here we would note that singularities appear in two places here. Of course, there is the 
geometric singularity in the entangling surface ci. However, the conformal transformation 
taking us from the flat metric (A.l) to the cylindrical metric (A. 2) is also singular at precisely 
the same point, i.e., this transformation maps the origin in R'^ toY ^ —00 in i? X 5 . 
Hence it seems this transformation is 'anomalous' in that it does not preserve the coefficient of 
the universal contribution to the entanglement entropy. This effect is somewhat reminiscent 
of the anomaly that arises in mapping a straight Wilson line in N=4 SYM to a circular Wilson 
loop, both in R^. While the former has a vanishing expectation value, the expectation of the 
latter yields a nontrivial result [43, 44]. 
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As a final example, we consider the EE for a cylinder in i? x S^. In this case, we insert 
d = 5 into eq. (A. 9) and follow the steps in section 3.1 to find that 



55 1 cylinder 



4 



L^sin^l^ L^sin^l^ 

log{H/S) + log(/io) 



3(53 
4 L cos^ Q 
95 



3<53 

log(/io) + <?5log(V^) + 0(5°) 



4 L cos^ $7 
95 



log(i7/5) 



(A.13) 



where precisely the same as given in eq. (3.26). Hence in comparing eqs. (A.13) and (3.25), 
we clearly see that the universal logarithmic term for a cylinder in i? x 5^ matches with that 
for cone C2- However, as compared to EE for C2 in (B.l), there are again new non-universal 
divergences, which take the form log(5)/53 and log(5)/5. 

We conclude this section by saying that although we have studied examples in d = 3, 4 
and 5, we expect that the log^ S term for cone Cd~3 in even d is not invariant under the 
conformal mapping from flat space to a cylinder. However, for odd d, the log 5 term will be 
invariant under the corresponding conformal transformation. 



B. Intermediate quantities for calculation of EE 

In this section, we mention various intermediatory steps in calculation of EE 
B.l EE for cone C2 in d = 5 

\7r^ L^lH^sin'^in) 4Hcos'^{n) 



S5 



C2 



+ q^\og{6/H) + 0{5'' 



(B.l) 



953 95 

where is given by (3.26). The h{6) in (3.26) is solution of following equation of motion 

h{l + /i^) sm{e)h + 2h cos{e)h^ + 2 (2 + h^) s\n{e)h^ 

+2h{l + h^) cos{0)h + (4 + Ih^ + 3h^) sm{e) = . (B.2 

B.2 EE for cone in d = 6 

TT^ I H'^sin^in) 27i?2cos2(J7)sin(0) 9 cos(17) cot(17)(31 - cos(2J7)) 



5, 



1C3 



1654 



12852 



+ 



8192 



/ sin3(17) 27cos2(17)sin(17) 



9cos(0) cot(fi)(31 -cos(20))log(/io) 
4096 



log(5/F)^ 
(B.3) 



log (5/i/) + 0(5° 



where q^ is given by 
Jo 



dh 



sin3(0) r ~ ~ sin2(fi) _ 27cos2(f^)sin(fi) 



+ 



9 cos(17)(31 - cos(20)) cot(O) 
4096/1 



(B.4) 
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B.3 EE for conical singularity ci x i?^ in d = 6 

For this case, the integrand in EE behaves as 

sm{0)-\/i+h^+h^ sin(f2) 3cos(il) cot(ri) 3(13 — 19cos(2il)) cot^(fi) csc(fi) 

hhP hp- + ^ 4096/1 ^ 

cot2(0) csc3(f^)(-375 + 348cos(2Q) -45cos(4fi) - I2863 sin(2f^) + 6463 sin(4J^)) 

32768 

+0{h), 

near the boundary. Using this, we find that EE is given by 



hlog{h) 
(B.5) 



Cl 



H^sm{n) 3cos(f])cot(0)log(5) 



8S^ 



+ 



64(52 



3^^(13- 19 cos(2J^))cot2(n) csc(O) + 384^gcos(J^) cot(O) (1 - 21og(/io)) - 4096sin(J^) 



16384/i2,52 



1 r° rln .2 T (n^ M 3(13 - 19 COs(2^])) COt' {Q) CSc(^^) Iog(<5) o^ 

■2^70 ^^^-^^(^^^ 8192^^ 



(B.6) 



where 

J6(/i) ^ 



'sin(^)Vl + ^2 + /j2 ^ gin(j^) ^ 3cos(Q)cot(n) ^ 3(13 - 19cos(2n)) cot2(Q) csc(Q)' 



hh^ 



32^3 



4096^ 



(B.7) 



B.4 EE for cone C2 in CFT dual to Gauss-Bonnet gravity 

For five-dimensional cone i?+ x in CFT dual to the Gauss-Bonnet gravity, the entanglement 
entropy turns out to be 



55 = 

^IC2 



H^sin^in) (1 - 4A/oo) 2Hcos'^{n) (2 - 7XU 



9(53 



9(5 



(B.8) 



(2cos2(J^) (2 - 7A/oo) hi - 3sin2(J^) (1 - 4A/oo)) 



9hl 



where 



£3 = 3sin^(e) (1 + /i^ + ^A}?f^(h^ (2 + S/i^ + h^) sin^(0) 

+4/t^ (2 + 3/i^ + /i^) cos(^) w^{ff)h + /i^ (4 + 6/1^ + 3/i^ + (8 + 9/1^ + 3/i^) cos(2e)) €\^{Q)h^ 
-h (4 - /i^ - (4 + 3/1^ + cos(20)) sin(20)/i3 + (/i^cos^(^) - 3h^ cos'^{e) sm^{e) + 6sin^(0)) /i^) 

-2Asin3(^)/oo (l + h^ + /i^) (-4/icos(e)/i + sin(0) (g + + 9^^)) 

£4 = 3^^AY^l + /i2 + /i2^sin2(e) (l + h^ + A^) + 2A/oo( (2 + 3/1^ + /i^) sin2(e) 

+/i (2 + h^) sm{2e)h + (/i^ cos2(0) - sin2(0)) /i^] \ . (B.9) 
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B.5 EE for cone C3 in CFT dual to Gauss-Bonnet gravity 

For cone i?"*" x S*^ in d = 6 dimensional CFT, which is dual to Gauss-Bonnet gravity, the EE 
is given by 

Att^L'' ri^'^sin3(r2)(3- lOA/00) if^ cos2(0) sin(f^) (27 - 86A/00) 



C3 



48^4 



128(52 



, 3 cos(n) cot(O) ((93 - I9OA/00) - (3 - 2A/00) cos(20)) , 2 

H t::^^:^ ^og[d/H) 



8192 



1 



12288/1^ 



- 27h^ cos(J^)(31 - cos(2^^)) cot(^^) log (ho) - 5184/ig cos^(fi) sin(J^) 
+3072 sin2(Q) -2A/00 (9/1^ cos(Q)(-95 + cos(20))cot(0) log (/lo) (B.IO) 



-8256/i§ cos2(^^) sin(i^) + 5120sin^(J)) + / dh{U- 4) \og{5/H) + C)((5°) 



ho 



Here i^s and are terms which come from making the h integrand in original expression of 
EE finite and these are given by 

£5 = ^ with (B.ll) 

JC5N = sin ^ ^12 A V^/i^ sin^ ^ + 48 A cos 6 h' sin^ ^ /t + 24 A/00 cos 6 sin^ ehh{l + {l- 5\f^) /i^) 

+6 A fooh^ sin(20) h (sin^ 9 + 2\f^(5 sin^ + 2 cos^ 9 h^^^ + \ f^h^ sin^ A sin^ (3 + A^) 

+2 A /oo (18 sin^ - (1 + 17 cos(2e))/i2) j + 2 A /oo/i^ (s sin^ O + Xf^ (l4 sin^ + 9 sin2(20) /i^) ) 

+2 sin^ (^20AV^/i^ + 3(1 + h^f - 8 A /oo(l + 3 /i^ + 2 A^)) + h'^ sin^ ^12 sin^ 0(1 + /i^) 

+4 AV^A^ (7 + 17 cos(20) - 9 cos^ eh^)-X /oo (24 sin^ + (17 - 23 cos(20))A^ + 6 cos^ 
+/i^ (^6 sin^ + A /oo sin^ (l + cos(20) + Ql?) + 4 A^^ (4 sin^ + 3 (5 + 8 cos(20)) sin^ f? 

+3 cos^0/t^) n , 



Ad = 6 /i^ /lA/l + ^^ + ^^^^/oo/i sin^^ + 9A/ooCOS0 sin 0/1/1 + 6A/00 cos sin 0/i^/i 



and 



+ sin^ (^1 + + 2 A /oo (2 - /i^)) + /i^ (sin^ + A/00 (7 sin^ + 3 cos^ }? 



sin^ Q (3 - lOA /oo) cos2 sin J] (27 - 86A/0 



3/i5 



+ 



32/i3 



3 cos O cot O (93 - 3 cos{2Vt) + 2 A (-95 + cos(20))/^ 



4096/1 

Note that in the asymptotic limit, h ^ Q and = C'^ + 0{h). 



(B.12) 
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B.6 EE for k x S"^ in CFT dual to Gauss-Bonnet gravity 

The values of £'s in eqn. (5.28) are following: 

Co = . . ^( (1-6/ooA + 8/^A2)A^ 

hWl + h^ + h^ ((1 _ 2/^A) /i2 + (1 - 2/ooA) + 4/ooA + 1 



N 



+2 (1 - 5/ooA + (1 - 6/ooA + 8/^A2) h^) 

+ (1 - 6/ooA + 8/^A2) h^ + {2- 10/ooA) - 4/ooA + l) , 

2/ooAVi + /i2 + A2 
3P 



72/l2 (1 + /i2) (l + /i2 + /,2^ 3/2 (^(1 _ /l2 + (1 - 2/ooA) 7^2 + 4/ooA + 1^ 

5 (3 - 10/ooA + 8/^A2) (1 + h^) + (42 - 26/ooA - 56/^A2 + 3 (29 - 80/ooA + 64/^A2) 

+15 (3 - 10/ooA + 8flX^) h^)h^ + (39 + 8/ooA + 48/^A2 + (l23 - 172/ooA + 64/^ A^) 

+3 (43 - 110/ooA + 88/^A2) + 15 (3 - 10/ooA + 8/^A2) /i^) + 5 (3 - 10/ooA + 8/^A2) 
+ (57 - 140/ooA + 112/^A2) /i^ + (81 - 146/ooA + 120/^A2) 

+ (51 - 72/ooA + 48/^A2) - 16/ooA + I2' 

Cb = z (-C2P + C3P - C3P + £5p) 

- 2/ooA) g^h^h + Q^h ( 1 + 6/ooA + (1 - 4/ooA) h^] 

(B.13) 



2/ooA/ig3 + 2/ooA/i3^3 + (1 - 2/ooA) g3h^h + g^h | 


(l + 6/ooA + (l-4/ooA)A2) 


1 


fl + /l2 + hA 


1 3/2 
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